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KM  =  16,  Gaussian  Signal .  32 

5  Required  Input  SNR  p(dB)  for  Pf  =  IE-6,  Pd  =  0.9, 
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6  Required  Input  SNR  p(dB)  for  Pf  =  IE-3,  Pd  =  0.5, 
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KM  =  256,  Gaussian  Signal .  33 
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H-l  Estimated  Covariance  Coefficients  for  K  =  4 .  H-2 
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DETECTION  PERFORMANCE  OF  OR-ING  DEVICE  WITH 
PRE-  AND  POST-AVERAGING:  PART  I  -  RANDOM  SIGNAL 

INTRODUCTION 

The  need  to  process  and  condense  large  amounts  of  data  is 
encountered  frequently  in  modern  Navy  systems  that  employ 
multiple  beams,  frequency  bins,  range  cells,  et  cetera.  One  way 
of  accomplishing  this  goal  is  by  or-ing  a  number  of  inputs  into  a 
single  output;  that  is,  allow  only  the  largest  of  a  set  of 
quantities  to  pass  on  for  further  processing  and  completely 
reject  the  remainder.  However,  since  this  or-ing  operation  is 
highly  nonlinear,  destroys  information,  and  tends  to  cause  small- 
signal  suppression,  some  pre-averaging  of  the  inputs  to  the 
or-ing  device  is  often  employed  in  an  effort  to  build  up  the 
signal-to-noise  ratio  (SNR)  prior  to  the  maximum  comparison. 
Additionally,  at  the  or-ing  output,  some  additional  post¬ 
averaging  is  frequently  employed,  once  again  in  an  effort  to 
build  up  the  SNR,  this  time  before  a  threshold  comparison  is  made 
for  purposes  of  declaring  a  signal  present  versus  absent.  The 
pertinent  block  diagram  is  displayed  in  figure  1. 

There  are  N  channels  of  real  input  data  available  for 
processing,  namely,  {xn(t)}  for  1  <  n  <  N,  where  time  has  been 
normalized  so  that  time  sampling  instant  t  is  integer.  Under 
hypothesis  HQ,  there  is  Gaussian  noise  only  in  all  the  inputs, 
whereas  under  hypothesis  H^,  a  signal  is  also  present  in  one 
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PRE-  POST¬ 
INPUTS  AVERAGER  OR-ING  AVERAGER  OUTPUT 

Figure  1.  Or-ing  with  Pre-  and  Post-Averaging 

(unknown)  channel.  The  goal  of  the  processor  in  figure  1  is  to 
determine  signal  presence  with  a  high  detection  probability  while 
realizing  a  specified  acceptable  low  false  alarm  probability. 

Each  pre-averager  accumulates  K  statistically  independent 
consecutive  samples  of  its  corresponding  input  xn(t),  yielding 
output  yn ( t ) ,  which  is  then  subjected  to  or-ing  amongst  N 
competitors.  The  or-ing  output  is 

v(t)  *  max  { y^  ( t  ),*•*,  yN  ( t ) }  .  d) 

Finally,  the  post-averager  accumulates  M  samples  of  its  input 
v(t),  and  compares  its  output  w(t)  with  a  fixed  threshold. 

There  is  an  additional  complication  of  the  post-averaging 
procedure,  which  is  not  indicated  explicitly  in  figure  1.  In 
particular,  the  or-ing  output  v(t)  is  sampled  every  J-th  time 
instant  and  summed  over  M  samples  to  give  processor  output 
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M-l 

w(t)  =  )  |  v(t-mJ)  =  v(t)  +  v(t-J)  +  •••  +  v(t-J(M-l)) 

m=0 

=  maxty^t)  yN  ( t ) }  +  maxty-^t-J)  ,...,  yN(t-J)}  +  ••• 

+  max{y1(t-J(M-l) )  yN( t-J(M-l ) ) }  ,  (2) 

where  J  is  the  skip  factor.  If  skip  factor  J  =  K,  then  system 
output  w(t)  is  the  sum  of  M  nonoverlapping  adjacent  blocks  of 
pre-averaged  data,  each  of  length  K,  covering  a  total  of  KM  time 
instants . 

More  generally,  for  arbitrary  J,  since  the  pre-averager  spans 
a  time  interval  of  length  K-l,  and  the  post-averager  spans  time 
interval  J(M-l),  the  total  observation  time  is  composed  of 

T  =  J(M-l)  +  K  (3) 

samples.  Alternatively,  the  last  term  in  equation  (2)  involves 
the  quantities  {xR( t-J(M-l )-K+l ) } ,  whereas  the  first  term  in 
equation  (2)  involves  {xn(t)},  for  a  total  of  T  time  instants. 

More  accurately,  T  can  be  interpreted  as  the  time-bandwidth 
product  of  the  input  data  utilized  by  the  processor  in  figure  1. 
Hypothesis  H1  presumes  that  a  signal  is  present  in  the  same 
(unknown)  channel  for  time  instants  t,  t-1,...,  t-J(M-l )-(K-l ) , 
for  a  total  of  T  samples.  If  J  =  K,  then  equation  (3)  reduces  to 
time-bandwidth  product  T  =  KM,  as  above;  this  special  case  is 
called  block  processing  and  is  the  major  focus  of  this  study. 
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The  reason  for  allowing  1  <  J  <  K  is  that  the  possibility 
exists  that  some  information  may  be  lost  unless  intermediate 
values  of  or-ing  output  v(t)  are  observed,  at  least  more 
frequently  than  every  K  samples;  see  equation  (2).  However,  this 
general  case  causes  statistically  dependent  non-Gaussian  random 
variables  (RVs)  to  be  accumulated  in  the  post-averager,  which  in 
turn  leads  to  an  intractable  analysis  problem.  Therefore,  the 
study  of  the  detectability  performance  for  this  particular 
feature,  namely,  skip  factor  J  <  K,  must  be  accomplished  by 
simulation.  Even  if  every  xn(t)  is  statistically  independent  of 
every  xm(u)  if  n  /  m  or  t  /  u,  the  system  output  w(t)  is  a  sum  of 
statistically  dependent  RVs  with  non-Gaussian  statistics,  when 
J  <  K.  Finally,  the  case  of  J  >  K  is  not  considered  at  all, 
because  that  leads  to  discarding  relevant  new  information  in 
or-ing  output  v(t)  and  a  definite  degradation  in  performance. 


SCALING  PROPERTY 

Consider  a  set  of  scaled  inputs  (xn(t)}  such  that 

xn(t)  =  a  xn(t)  +  b  ,  a  >  0  ,  1  <  n  <  N  .  (4) 

It  can  then  be  easily  shown  that  the  corresponding  pre-averager 
outputs  Yn(t)  =  a  Yn(t)  +  bK,  and  that  the  or-ing  outputs  satisfy 
v(t)  «*  a  v(t)  +  bK,  leading  to  output  relationship  w(t)  =  a  w(t) 

+  KM.  Thus,  w(t)  is  a  scaled  version  of  w(t),  along  with  an 
additive  constant.  Since  the  transformation  between  the  two 
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processor  outputs  w(t)  and  w(t)  is  monotonic,  the  receiver 
operating  characteristics  (ROCs)  of  the  two  or-ing  processors, 
one  operating  on  {xn(t)}  and  the  other  operating  on  {xn(t)},  are 
identical,  regardless  of  the  values  of  scaling  parameters  a  and  b 
in  equation  (4).  Then,  without  loss  of  generality,  the  scalings 
can  be  chosen,  for  example,  so  that  input  xn(t)  has  zero  mean  and 
unit  variance  under  Hq,  if  desired.  That  still  allows  xR(t)  to 
have  arbitrary  mean  and  variance  under  H^.  Alternatively,  if 
xn(t)  is  always  a  positive  RV,  the  choice  of  b  =  0  with  unit 
variance  of  xR(t)  under  Hq  is  a  possibility.  This  still  leaves 
the  variance  of  xn(t)  under  arbitrary.  The  point  is  to  choose 
the  two  free  parameters  a  and  b  in  equation  (4)  so  that  the 
statistics  of  {xn(t)}  under  Hq  are  as  simple  as  possible;  this 
feature  will  be  utilized  without  comment  in  the  future 
developments . 


STANDARD  OF  COMPARISON 

The  standard  of  comparison  against  which  to  measure  system 
performance  losses  is  N  =  1,  no  or-ing.  Then,  for  arbitrary 
pre-averager  weights  { g ^ }  and  arbitrary  post-averager  weights 
{ hm> ,  the  various  processor  waveforms  for  N  =  1  in  figure  1  are 

y-l^t)  =  Y  9k  xi ( t-k )  ,  V(t)  -  yx(t)  , 

k 

w(t)  =  Y  hm  v(t-m)  =  Y  cj  /  (5) 

El  j  J 
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where  weights  {c^}  are  the  convolution  of  the  pre-averager 
weights  { g ^ }  and  the  post-averager  weights  {hmJ  (including  skip 
factor  J).  For  statistically  independent  time  inputs  {x^(t)}, 
the  characteristic  function  (CF)  of  processor  output  w(t)  is 
given,  for  N  =  1,  no  or-ing,  by  ensemble  average 

fw(0  =  exp(Uw(t))  =  TT  fv(c,U  ,  (6) 

W  j  A  J 

using  equation  (5),  where  fx  is  the  CF  of  x^(t)  under  the 
appropriate  hypothesis,  either  Hq  or  .  This  particular 
processor  output  CF  (for  N  =  1)  can  be  accomplished  analytically 
and  then  Fourier  transformed  numerically  to  get  accurate  ROCs  for 
any  K,  M,  and  sets  of  weights  {g^}  and  { hm } . 


OUTLINE 

For  skip  factor  J  =  K,  namely,  block  processing  in  figure  1, 
expressions  will  be  derived  that  allow  for  accurate  evaluation  of 
the  false  alarm  probability  and  the  detection  probability 
for  the  decision  variable  w(t).  Furthermore,  this  will  be 
accomplished  for  arbitrary  amounts  of  pre-averaging  K,  arbitrary 
amounts  of  or-ing  N,  arbitrary  amounts  of  post-averaging  M,  and 
arbitrary  input  SNRs.  No  approximations  are  involved,  the 
analysis  is  not  limited  to  mean  and  variance  calculations,  and  no 
appeal  is  made  to  the  central  limit  theorem.  Rather,  the 
approach  employs  a  judicious  combination  of  analysis  with 
computer-aided  numerical  calculations.  The  accuracy  of  the  end 
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result  is  limited  only  by  the  computer  accuracy.  An  entire  ROC 
can  be  generated  in  minutes. 

There  are  three  different  input  signal  models  in  additive 
Gaussian  noise  that  are  of  interest  here;  they  are: 

(I)  Random  (Gaussian)  signal, 

(II)  Phase-incoherent  signal,  and 

(III)  Coherent  (deterministic)  signal. 

This  report  will  address  all  three  signal  models,  but  will 
concentrate  numerically  on  part  I,  the  random  Gaussian  signal, 
and  will  present  numerous  ROCs  that  completely  characterize  the 
performance  of  the  or-ing  system  in  figure  1.  Some  related  past 
work  on  or-ing  is  listed  in  references  1  through  9. 

For  all  three  signal  models,  the  optimum  processor  can  also 
be  determined  and  its  ROCs  calculated.  Then,  the  exact  losses 
that  the  or-ing  procedure  causes  can  be  accurately  quantified. 
That  chore  will  be  completed  here  in  part  I  for  the  random 
signal;  parts  II  and  III  are  reserved  for  follow-on  work. 

The  analytical  approach  utilized  here  can  be  extended  to 
include  quantizers  on  the  inputs  { x  ( t ) }  in  figure  1  and  still 
obtain  accurate  ROCs  for  decision  variable  w(t).  In  fact,  not 
only  can  there  be  arbitrary  amounts  of  quantization  (L  levels), 
but  the  input  probability  density  functions  (PDFs)  can  be 
completely  arbitrary.  Documentation  of  this  work  is  reserved  for 
the  future. 
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Some  or-ing  processors  perform  quantization  on  the  output 
v(t)  of  the  Greatest-Of  device  in  figure  1,  prior  to  post¬ 
averaging.  The  PDFs  at  the  output  of  this  quantizer  are  then 
collections  of  impulses,  whose  areas  are  governed  by  the 
cumulative  distribution  functions  of  v(t)  under  hypotheses  Hq  and 
H^.  If  the  quantizer  levels  are  equally  spaced,  then  the  impulse 
locations  are  also  equally  spaced,  and  a  complete  analytic 
approach  is  possible,  all  the  way  to  processor  output  w(t)  in 
figure  1.  See  references  3  and  4,  for  example. 
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INPUT  DATA  MODELS 


In  the  following,  gk  and  hk  are  independent  zero-mean  unit- 
variance  Gaussian  RVs,  and  ek  is  an  exponential  RV  with  unit 
mean,  that  is,  with  PDF  exp(-u)  for  u  >  0. 


RANDOM  GAUSSIAN  SIGNAL  IN  GAUSSIAN  NOISE 


In  the  signal  channel,  the  system  input  to  figure  1,  at  time 
tk  under  hypothesis  H^,  is  an  envelope-squared  variate 


w 


-*[c. 


+  9u)  +  ( b. 


V 


(7) 


where  signal  components  ak  and  bk  are  independent  zero-mean 
Gaussian  RVs  with  common  variance  pk,  and  gk  and  hk  are  additive 
Gaussian  noise  components.  The  input  mean  is  then 


xk  =  1  +  pk  (8) 

for  signal  present.  Since  xk  is  an  envelope-squared  quantity, 

Pk  is  an  input  (per  sample)  power  SNR  measure. 


The  CF  of  xk  under  H1  is,  from  equations  (7)  and  (8), 

£k(»  -  e»p(Uxk)  -  mTil  T  >k)  •  (9) 

Since  this  CF  corresponds  to  an  exponential  RV  with  mean  1  +  pk, 
it  is  possible  to  use,  instead  of  form  (7),  input  data  value 
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(10) 


xk  =  U  +  Pk)  ek  , 

both  for  analysis  as  well  as  simulation  purposes  under  H^.  For 
the  noise  channels,  pk  =  0  and  xk  =  ek  under  HQ. 

When  pk  =  p  for  all  k,  that  is,  constant  signal  powers  with 
time,  the  pre-averager  output  for  the  signal  channel  under  is 

K  K 

y  =  )  1  x,  =  (1  +  p)  >  I  e.  »  (1  +  p)  y°  ,  (11) 

k=l  K  k=l 

where  y°  is  the  corresponding  noise-only  pre-averager  output. 

This  property  can  be  used  advantageously  for  both  analysis  as 
well  as  simulation  purposes.  The  ROCs  can  be  parameterized  by 
input  power  SNR  p;  this  is  the  mathematical  setup  for  part  I. 


ALTERNATIVE  ADDITION  OF  ENVELOPES  FOR  RANDOM  GAUSSIAN  SIGNAL 

If  the  input  data  value  in  figure  1  is  an  envelope  instead, 
denote  it  by  xk  =  (Xj^)*5  in  terms  of  the  RV  in  equation  (7).  This 
means  that  xk  can  be  generated  alternatively  according  to 

xk  =  (1  +  Pjj)^  •  (12) 

The  mean-square  value  of  envelope  xk  is  then 

x£  =  1  +  pk  ,  (13) 

where  pk  is  an  input  power  SNR  measure. 
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For  the  noise  channels,  pk  =  0,  thereby  yielding  envelope 


where  ek  is  a  unit-mean  exponential  RV. 


(14) 


When  pk  =  p  for  all  k,  the  pre-averager  output  for  the  signal 
channel  under  is 

K  K 

X  =  XI]  *k  -  U  +  p)H  XI]  ek  =  (1  +  p)H  1°  ,  (15) 

where  y°  Is  the  corresponding  noise-only  pre-averager  output. 

The  ROCs  can  be  parameterized  in  terms  of  input  power  SNR  measure 
p.  The  investigation  in  this  report  on  the  random  signal,  part 
I,  will  not  address  the  sum  of  envelopes,  as  given  by  equation 
(15),  but  will  concentrate  on  the  sum  of  envelopes  squared,  as 
given  by  equation  (11),  with  constant  powers  pk  «=  p. 


PHASE-INCOHERENT  SIGNAL  IN  GAUSSIAN  NOISE 

In  the  signal  channel,  the  input  for  this  signal  model,  part 
II,  at  time  tk  under  is  an  envelope-squared  variate 

xk  "  l[(9k  +  Ak  cos6k)2  +  (hk  +  Ak  ®inek)2]  - 

where  signal  amplitude  Afc  is  nonrandom,  and  phase  9k  is  a  RV  with 
arbitrary  PDF,  that  is,  incoherent  phase.  The  input  mean  is 
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(17) 


*k  -  1  +  \  Ak  E  1  +  pk  ■ 

regardless  of  the  value  of  ek-  Since  xk  is  an  envelope-squared 

2 

quantity,  pk  =  Ak/2  is  an  input  (per  sample)  power  SNR  measure. 

The  CF  of  xk  under  ,  conditioned  on  a  particular  value  of 
RV  0k,  is,  by  use  of  equation  (16), 

£kU|ek>  -  exp(Hxk)  ■=  exP (ek  1  •  (18) 

But,  since  the  right-hand  side  of  CF  (18)  is  independent  of  0k, 
it  is  possible  to  use  the  data  value 

xk  -  l[<*k  +  V2  *  hk]  <19) 

for  the  signal  channel  input  under  H^.  On  the  other  hand,  for 
the  noise  channels,  Ak  =  0  and  CF  (18)  reduces  to  (1  -  i£) 
which  corresponds  to  a  unit-mean  exponential  RV.  Therefore,  for 
the  noise  channel  inputs,  use  the  quantity 

xk  -  ek  •  (201 


When  Ak  =  A  for  all  k,  constant  signal  amplitudes  with  time, 
the  pre-averager  output  for  the  signal  channel  can  be  expressed 


K  .  K  -  ~  K  .2 

Y  =  II]  xk  =  J  £  ( gk  +  hk )  +  A  C  9k  +  K  f 

k— 1  k — 1  k — 1 

K  2 

■  y  +  A  c  gv  +  K  T  ,  (2D 

k=l  K  z 
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under  H^,  where  y°  is  the  corresponding  noise-only  pre-averager 

output.  (However,  the  k-th  term  in  the  sum  for  y  cannot  be 

2 

replaced  by  +  A  g^  +  A  /2,  where  is  an  independent 
exponential  RV,  because  the  CF  of  this  latter  combination  is 

1  A2 

— Z-Jl  explp  U(1  +  i^)]  '  P  *  §  /  (22) 

which  is  not  identical  in  form  to  equation  (18).  The  two  CFs 
do,  however,  have  the  same  mean  and  variance.)  The  ROCs  for 
figure  1  using  pre-averager  (21)  can  be  parameterized  by  input 
amplitude  SNR  measure  A  or  by  input  power  SNR  measure  p. 

Because  of  the  way  parameters  {p^}  appear  in  equation  (18), 
the  CF  of  the  pre-averager  output  can  be  found  in  the  closed  form 

fy(t>  =  u-urK  e*p(„K  jifj)  ,  pK  -  E  pk  • 

Therefore,  all  the  following  results  pertaining  to  the  phase- 
incoherent  signal  with  equal  SNRs  p  and  sums  of  squared-envelopes 
actually  hold  for  arbitrary  SNRs  {oy,} ,  if  p  is  replaced  by  pR/K. 
However,  this  observation  is  of  limited  utility  unless  all  of  the 
M  sums  of  {p^},  one  for  each  post-averaging  interval,  are 
identical . 

ALTERNATIVE  ADDITION  OF  ENVELOPES  FOR  PHASE- INCOHERENT  SIGNAL 

If  the  input  data  value  in  figure  1  is  an  envelope  instead, 
denote  it  by  =  (x^*5  in  terms  of  the  RV  in  equation  (16). 
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Since  the  conditional  CF  of  xk  is  independent  of  ©k,  according  to 
equation  (18),  the  conditional  PDF  of  xk  is  independent  of  ©k. 
Therefore,  any  nonlinear  transformation  of  x^,  such  as  x^,  will 
also  have  a  conditional  PDF  independent  of  ©k.  This  means  that 
xk  can  be  generated  alternatively  according  to 


C2  [( gk  +  Ak 


(23) 


under  H^.  The  mean-square  value  of  envelope  xk  is  then 

^k  “  1  +  I  Ak  s  1  +  pk  '  (24) 

where  pk  is  an  input  power  SNR  measure. 


For  the  noise  channels,  Ak  =  0,  thereby  yielding  envelope 


— k  ~ 


(25) 


where  ek  is  a  unit-mean  exponential  RV. 


When  Ak  =  A  for  all  k,  the  pre-averager  output  for  the  signal 
channel  under  is 

Y  -  C  2k  -  C  (|(9 i  +  hj)  +  A  9k  +  ^  .  (26) 

k — 1  k*l  ' 

which  cannot  be  simplified  or  expressed  in  terms  of  alternative 
RVs .  The  ROCs  can  be  parameterized  in  terms  of  input  amplitude 
SNR  measure  A  or  in  terms  of  input  power  SNR  measure  p  =  A2/2. 
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COHERENT  (DETERMINISTIC)  SIGNAL  IN  GAUSSIAN  NOISE 


In  the  signal  channel,  the  input  for  this  signal  model,  part 
III,  at  time  t^  under  is  amplitude  variate 


xk  =  *k 


+  d 


k  ' 


(27) 


where  d^  is  nonrandom.  The  input  mean  is  x^  =  d^,  while  the 

input  variance  is  var(x^)  =  1.  Thus,  d^  is  an  input  amplitude 

2  2 

SNR  measure,  while  dk  (or  d^/ 2)  is  an  input  power  SNR  measure. 

For  the  noise  channels,  x^  =  gk,  with  zero  mean  and  unit 
variance.  Thus,  dk  can  be  interpreted  as  a  per-sample  difference 
in  mean  inputs,  divided  by  the  common  input  standard  deviation. 


When  dk  =  d  for  all  k,  constant  signal  amplitudes  with  time, 
the  pre-averager  output  under  for  the  signal  channel  is 

K  K 

y  =  HU  xk  =  HH  gk  +  K  d  =  y°  +  K  d  ,  (28) 

k=l  k=l 

where  y°  is  the  corresponding  noise-only  pre-averager  output. 

The  ROCs  can  be  parameterized  by  input  amplitude  SNR  measure  d  or 
input  power  SNR  measure  d  (or  d  /2).  More  generally,  the  single 
parameter  value  dT  =  Ek  dk  characterizes  the  performance  of  the 
processor  when  the  amplitude  SNRs  {d^}  are  unequal  and  arbitrary. 
However,  this  observation  is  of  limited  utility  unless  all  of  the 
M  sums  of  {d^},  one  for  each  post-averaging  interval,  are 
identical . 
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OPTIMUM  PROCESSING  OF  AVAILABLE  DATA 


The  optimum  processors  for  all  three  signal  models  are 
derived  in  appendix  A  for  the  case  where  channel  occupancy 
identification  is  required,  as  well  as  the  alternative  case  where 
channel  identification  is  irrelevant.  For  all  situations,  the 
optimum  processor  is  essentially  given  by  threshold  comparison 


max 
l<n<N  U=1 


(n 


xn(t) 


> 

< 


(29) 


see  equations  (A-13),  (A-19),  and  (A-26).  This  operation  is 
equivalent  to  doing  all  pre— averaging  and  no  post— averaging  in 
figure  1,  that  is,  take  K  =  T  and  M  =  1.  Doing  so  defers  the 
nonlinearity  (or-ing)  until  after  all  possible  linear  combining 
(pre-averaging)  has  been  accomplished.  This  case  will  be 
thoroughly  numerically  evaluated  and  will  serve  as  the  basis  of 
comparison  for  the  various  combinations  of  pre-averaging  and 
post— averaging,  that  is,  H  >  1  or  K  <  T. 
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PERFORMANCE  ANALYSIS  FOR  BLOCK  PROCESSING 

Due  to  the  pre-averaging,  each  RV  yR(t)  in  figure  1  is  a  sum 
of  K  independent  RVs  with  identical  statistics.  Also,  all  N 
channel  inputs  are  statistically  independent  of  each  other. 

In  order  for  or-ing  to  be  present,  N  >  2  is  required  in  the 
following  considerations.  (The  special  case  of  N  =  1,  no  or-ing, 
has  already  been  addressed  in  equations  (5)  and  (6).) 

For  a  noise-only  channel,  let  Pq  be  the  PDF  of  pre-averager 
output  yn(t),  and  let  co  be  the  corresponding  cumulative 
distribution  function  (CDF).  For  signal  present  in  channel  1, 
let  pj^  and  c1  be  the  PDF  and  CDF  of  y^t),  respectively. 

The  or-ing  output  v(t)  in  figure  1  is 

v ( t )  =  max{y1( t) , . . . ,yN( t) }  .  (30) 

Its  CDF  for  signal  present,  that  is,  hypothesis  ,  is 

cyl(u)  =  c^u)  [c0(u)]N_1  ,  (31) 

which  leads  to  the  corresponding  PDF  of  v(t)  under  as 
pvl(u)  =  c^(u)  =  px(u)  [c0(u)]N_1  +  cx(u)  ( N— 1 )  lc0(u)]N“2  Pq ( u ) 
=  [c0(u)]N~2  [p1(u)  c0(u)  +  ( N— 1 )  p0(u)  c1(u)]  .  (32) 

The  signal-absent  PDF  of  v(t)  under  Hq  follows  immediately  as 
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(33) 


Pv0(u)  =  N  pQ(u)  [cQ(u) ]N  1  . 

The  CF  of  or-ing  output  v(t)  is  given  by  Fourier  transform 

fy(£)  =  exp (  i  £v  ( t ) )  =  J  du  exp(i^u)  py(u)  ,  (34) 

where  either  relevant  form,  p^Q  or  pyl  above,  is  to  be  used. 
Finally,  for  block  processing,  J  =  K,  there  is  no  overlap  of  the 
data  used  in  the  post-averager,  meaning  that  processor  output 
w(t)  is  the  sum  of  M  independent,  identically  distributed  RVs. 
Therefore,  the  CF  of  the  decision  variable  w(t)  is  simply 

fwU)  =  [fyU)]M  •  (35) 

One  final  Fourier  transform  is  required  to  manipulate  CF 
fw(£)  into  the  desired  exceedance  distribution  function  (EDF)  of 
processor  output  w(t);  see  reference  9.  This  EDF  will  be  the 
false  alarm  probability  Pf  or  the  detection  probability  Pd, 
depending  on  whether  PDF  py£)  or  pyl  is  used,  respectively. 

The  fundamental  Fourier  transform  in  equation  (34)  will 
generally  have  to  be  accomplished  numerically  by  means  of  a  fast 
Fourier  transform  ( FFT) .  If  pre-averager  output  statistics  pQ 
and  c0  can  be  found  in  closed  form  or  readily  computed  form,  this 
cascaded  FFT  approach  will  lead  to  exact  false  alarm  probability 
results  for  the  or-ing  processor  in  figure  1.  If  pre-averager 
output  probability  functions  p.^  and  c^  can  also  be  readily 
evaluated,  accurate  detection  probability  results  can  be 
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calculated  as  well.  Some  important  numerical  considerations  for 
this  cascaded  FFT  approach  are  presented  in  appendix  B. 

For  the  special  case  of  M  «  1,  no  post-averaging,  output  RV 
w(t)  =  v(t),  and  the  corresponding  EDFs  follow,  for  any  K,  N,  as 

ewi<”)  =  1  -  c1(u)  tc0(u)]N-1  ,  (36) 

ew0(u)  "  1  "  tco(u)]N  ‘  (37) 

Now,  return  to  arbitrary  values  of  K,  N,  and  M,  and 
specialize  the  general  results  above  to  the  three  signal  models 
of  interest.  J  is  set  equal  to  K,  that  is,  block  processing. 


COHERENT  (DETERMINISTIC)  SIGNAL 


Define  d  =  /K  d  and  auxiliary  Gaussian  functions 


<j>(x)  =  (2n)  35  exp(-x^/2)  ,  *(x)=Jdu$(u). 


(38) 


Then,  for  constant  signal  amplitudes,  the  requisite  probability 
functions  for  pre-averager  outputs  {yn(t)}  for  any  K,  all  u,  are 


P0(u) 

Px(u) 


1 

(2nK)Jj 

1 

(2nK  )J$ 


(39) 


(40) 
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c0(u)  -  *(t!) 


C1<“»  -  *(tI-  a)  • 


(41) 


/ 


Using  equations  (33)  and  (32),  the  PDFs  for  orbing  output 
v(t)  for  any  K  and  N,  under  Hq  and  H^,  are,  respectively, 

pv0(u)  =  A  *(tk)  *(tk)  '  (42) 

pvl(u)  =  A  * (tk)  KtI  -  *(tk)  +  (N_1)  *[A)  $(tk  -  ^)]  • 

(43) 

The  corresponding  CFs  are 

fv0(O  -  N  J  dx  exp( i^K^x)  *(x)  $(x)N_1  ,  (44) 

fvlU)  =  Jdx  exP<i^KJ5x)  *U)N~2  (+(x-d)  i(x)  +  ( N-l )  +  (x)  $ ( x-d )  ] . 

(45) 

When  M  =  1,  the  EDFs  in  equations  (36)  and  (37)  can  be  combined 
with  the  results  in  equation  (41)  to  immediately  yield  closed 
form  expressions  for  and  P^  for  any  K,  N,  and  SNR  p. 

For  the  special  case  of  N  =  2,  these  integrals  can  be 
evaluated  in  terms  of  the  error  function  of  complex  argument, 
w(  );  see  reference  10,  chapter  7.  There  follows 

fv0K)  =  exp(-r2)  w(y)  ,  (46) 

and 

fvi<*>  *  \  exp(-Y2-v2+i2Yv)  tw(Y+iv)  +  w(Y-iv)]  ,  (47) 
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k  k 

where  y  =  K^£/2,  v  =  K^d/2 .  The  amount  K  of  pre-ave raging  is 
arbitrary.  The  CFs  of  decision  variable  w(t)  are  given  by  the 
M-th  powers  of  the  above  two  results  under  Hq  and  H^, 
respectively,  when  N  =  2. 

For  N  =  1 ,  no  or-ing,  and  K  and  M  arbitrary,  these  results 
simplify  to  system  output  PDF 

pwl(u)  =  - ^~k  exp("  (U  ^TT-d)2)  »  T  =  K  M  ,  (48) 

wi  ( 2nT)^  v  ' 

and  output  EDF 
00 

pd(v)  =  J  du  pwl(u)  =  $(qd  -  ,  q  =  t*5  =  (KM)*5  .  (49) 

v 

Obviously,  for  N  =  1 ,  only  the  product  KM  matters,  at  least  for 
unity  weights  in  both  the  pre-  and  post-averagers. 


PHASE- INCOHERENT  SIGNAL 


Some  of  the  following  was  presented  in  reference  3,  pages 
52  and  53.  Define  £  =  K  p  with  (reference  10,  equation  6.5.11) 


en(x) 


(50) 


Then,  for  constant  signal  amplitudes  and  u  >  0,  the  pertinent 
PDFs  and  CDFs  of  { y n ( t ) }  are 
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P0<u)  =  uK  |k-1)!  ~  '  co(u)  =  1  “  exp(_u)  eK-i(u)  '  (51) 


K-l 

p^fu)  =  (j)  2  exp(-u-£)  xk_i  ,  (52) 

-  (1  -  U)"K  exp(p  Y^n)  r  (53) 

c1(u)  *  1  -  Qk((2£)Js,  (2u)Jsj  .  (54) 

The  PDFs  for  or-ing  output  v(t)  for  u  >  0  are  then 

pv0(u)  =  N  ^ — (K!lfi~U)  [1  "  exP(-u)  eK-l(u)lN_1  '  (55) 

K-l 


Pvl(u)  -  [1  -  exp(-u)  eK_1(u)]N  2  [(*-)  2  exp(-u-£)  xK_i  (2 (£u) **) 

x  [1  -  exp(-u)  eK_1(u)]  +  ( N-l )- — j  ^  { —  (l-QR  ( (  2£ )  ^  f  ( 2u ) %]  )  ]  . 

(56) 


in  the  special  case  of  R  *  1,  that  is,  no  pre-averaging,  the 
CF  of  v(t)  with  or-ing  and  signal  present  is  available  in  closed 
form  by  Fourier  transforming  pv^(u)  to  obtain,  using  reference 
11,  equation  6.631  4,  and  reference  12,  equation  (24), 


fvl<*> 


(-1) 


k+1 


U)(k  +  1  -  i£) 


exP("p  k  +  1  -\J(57) 


for  any  N  >  1.  The  mean  of  v(t)  with  signal  present  follows  as 
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g  (V)  kfrrir  «p(-  dhr) 


(58) 


Finally,  when  p  =  0,  equation  (57)  reduces  to  the  more  useful 
form  for  the  CF  as 


r  N  (  .  ..,-1  N  . 

£vo<«>  -  M1  -  *9]  >  'vo  ■  C  H  .• 

n=i  n=i 


(59) 


this  alternative  CF  form  will  be  derived  in  the  next  subsection. 


For  the  special  case  of  K  =  2  and  no  signal  present,  the  CF 
of  v(t)  is  available  by  Fourier  transforming  PDF  Pvq(u)  to  get 


N-l  ,  ,,n  n  ,  , . 

f  m  =  n i  \  1  (~1) \ — 1  _ *±1 _ 

v0vt*'  w*  — 1  (N-l— n)  i  L — >  k+2 

vu  n=0  1  n;*  k=0  (n-k)I  (n+l-H) 


(60) 


This  result  holds  for  arbitrary  amounts  of  or-ing,  N  >  1 


For  N  =  2,  the  CF  of  v(t)  under  is  available  (after 
considerable  labor)  by  Fourier  transforming  PDF  pv^(u)  to  get 


z,v  K-l 


$  if)  -  if)  £  W  -i  <“> 


for  any  K  >  1,  where  z1  =  1  -  i£,  z2  =  2  -  i£.  For  K  -  2,  the 
mean  follows  as  =  2(1  +  p)  +  exp(-p)  (6  +  p)/8,  using  £  =  2p. 
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For  N  =  2  and  signal  absent,  the  Fourier  transform  of  PDF 


Pv0(u)  yields  CF 


K-l 


£v0<<> 


(1  -  it)' 


-  2 


c  (K  -  £ +  k) 

k=0  v  K  J 


(2  -  H) 


K+k 


(62) 


with  mean 


K-l 


"vo  - »  (i  -  c  p : k)  ^rr) . 


(63) 


This  form  for  CF  fvg(ii)  in  equation  (62)  is  not  obtained  by 
setting  £  =  0  in  equation  (61);  however,  numerical  investigation 
of  both  forms  reveals  they  are  equal.  For  K  =  2,  the  CF  in 
equation  (62)  further  specializes  to 


fvO<S> 


2 _ 2 _ 4 

(1  -  U)2  (2  -  U)2  (2  -  U)3 


(64) 


The  CF  of  decision  variable  w(t)  for  signal  present  is  given 
by 

f„i(5>  -  t£vl<tn"  ■  (65) 

When  CF  fw^(S)  is  available,  a  Fourier  transform  will  yield  the 
detection  (and  false  alarm)  probabilities  for  any  K,  N,  and  M. 
When  M  =  1,  the  EDFs  in  equations  (36)  and  (37)  can  be  combined 
with  results  (51)  and  (54)  to  immediately  yield  closed  form 
expressions  for  Pf  and  Pd  for  any  K,  N,  and  SNR  p. 
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When  N  =  1 ,  no  or-ing  and  K  and  M  are  arbitrary,  the  results 


simplify,  respectively,  to  PDF  and  EDF 


T-l 

Pwl<u)  =  (|j  2  exp(-u-B)  IT_X  ^(Bu)*5)  for  u  >  0  ,  (66) 

00 

pd(V)  =  J  du  Pwl(u)  =  Qt  ((2b)15,  (2v)Jsj  for  v  >  0  ,  (67) 

v 


where  T  =  KM,  B  =  Tp  =  KMp.  Obviously, 
product  KM  matters;  this  presumes  unity 
and  post-averagers. 


for  N  =  1  here,  only  the 
weights  in  both  the  pre- 


RANDOM  GAUSSIAN  SIGNAL 

Some  of  this  material  appeared  in  reference  3,  pages  51  and 
52.  Define  a  =  1/(1  +  p).  Then,  for  constant  signal  powers  and 
u  >  0,  the  pertinent  PDFs  are 

K-l  . 

p0(u)  "  U  (K-l)!  '  c0(u)  "  1  -  exp(-u)  e^fu)  ,  (68) 

K  K—  1 

Pj^u)  -  - — U  (K-1^(~^U)  '  ^(u)  -  1  -  exp(-au)  eR_1(au)  .  (69) 

The  PDFs  for  or-ing  output  v(t)  are,  under  hypotheses  Hq  and 
H^,  respectively, 

K— 1 

Pv0(u)  =  N  U  (K!lH~U)  11  -  exP(-«)  eK-l(u)]N_1  '  <70> 

which  behaves  as  uKN-^  as  u  ->  0+,  and 
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X 


Pvl(u)  =  [1  -  exp(-u)  eK_1(u)]N  2  (K-l ) !  (aK  exP(“au) 

[1  -  exp(-u)  eK_i<u)l  +  (N-l)  exp(-u)  [1  -  exp(-au)  eR_1(au)]j. 


(71) 


In  the  special  case  of  K  =  1,  that  is,  no  pre-averaging,  the 
CF  of  v(t)  with  or-ing  and  signal  present  is  available  in  closed 
form  by  Fourier  transforming  PDF  pv^  to  get  (reference  11, 
equations  3.191  3  and  8.384  1) 

fvlU)  =  (1)N-1  ((a-U)N  +  (l-i^J  '  (72 

for  any  N  >  1.  Here,  (z)n  =  z(z+l)  •••  (z+n-1)  is  Pochhammerfs 
symbol;  see  reference  10,  equation  6.1.22.  When  p  =  0,  then 
a  =  1  and  this  CF  reduces  to 


(1) 


N 


(1-U) 


N 


(73) 


The  CF  of  decision  variable  w(t)  for  signal  present  is  given 
by  closed  form 

*wlU>  -  [«vl(?)]M  (74) 

in  this  special  case  of  K  =  1,  where  CF  fv^  is  given  by  equation 
(72).  Then,  one  single  Fourier  transform  will  yield  the 
detection  and  false  alarm  probabilities  (P^  and  P^)  for  any  N  and 
M,  that  is,  arbitrary  amounts  of  or-ing  and  post-averaging. 
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When  M  «  1,  the  EDFs  in  equations  (36)  and  (37)  can  be  combined 
with  results  (68)  and  (69)  to  immediately  yield  closed  form 
expressions  for  and  for  any  K,  N,  and  SNR  p. 


The  cumulants  {^(k)}  of  v(t),  for  signal  present,  are  rather 
involved,  even  for  K  =  1.  The  mean  and  variance  are  presented 
below,  while  the  moments  and  cumulants  of  general  order  k  are 
derived  in  appendix  C.  For  K  -  1, 


V1* 


N-l  , 

y  I 

n=l  n 


(1) 


+  r 


with 


r  = 


N-l 


(a) 


N 


(75) 


V2)  'Cl  -  (r  -  |)2  -  O  n(n+a-l )  •  (76) 

n=i  n  n=i 

For  signal  strength  p  =  0,  a  »  1  and  the  k-th  cumulant  of 
v(t)  is,  for  K  -  1, 

N  i 

Xv(k)  =  ( k-1 ) !  y  -f  for  k  >  1  .  (77) 

v  n=l  nK 

These  results  allow  for  arbitrary  amounts  of  or-ing,  that  is, 

N  >  1. 


For  the  special  case  of  K  =  2  and  no  signal  present,  the  CF 
of  or-ing  output  v(t)  is  available  in  closed  form  by  Fourier 
transforming  PDF  pv^(u)  to  get 


fv0<*> 


N-l 

=  Ni  y 

n=0 


(-1) 


n 


(N-l-n) 


n 

rz 

k=0 


k+1 


(n-k)!  ( n+l-i^ ) 


k+2 


(78) 
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This  result  holds  for  arbitrary  amounts  of  or-ing.  CF  fv^(£)  for 
K  =  2  can  also  be  found  in  closed  form,  although  it  is  more 
complicated  than  the  result  above.  Larger  values  of  K  yield 
progressively  more  complicated  closed  forms  for  fvg(S)  and 
fviU).  For  example,  pv^(u)  f°r  K  =  2  and  N  =  2  can  be  Fourier 
transformed  to  yield  CF 


fvlU> 


with  mean 


1  a2  (l+a2)U  -  (1+a)3 

2  +  2  +  3  ' 

(1  -  U)  (a  -  ur  (1  +  a  -  ur 


(79) 


Vl 


=  2 


2  3  4 

1  +  3a  +  3a  +  3a  +  a 

a  (1  +  a)3 


=  2 


11  +  22  p  +  18  p2  +  7  p3  + 

(2  +  p)3 


(80) 


When  N  =  1,  no  or-ing  and  K  and  M  are  arbitrary,  these 
results  simplify  to  processor  output  PDF 


wl 


(u>  - aT  “T^_nr-au)  *« «  >  ° 


(81) 


and  edf 

QO 

pd( v)  “  J  du  PW1^U^  “  exp(-av)  e^j^tav)  for  v  >  0  ,  (82) 

v 

where  T  =  KM.  For  N  =  1,  only  the  product  KM  matters;  the  above 
results  presume  unity  weights  in  both  the  pre-  and  post¬ 
averagers  . 
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QUANTITATIVE  PERFORMANCE  RESULTS 


This  section  gives  specific  quantitative  detectability 
results  for  or-ing  with  various  amounts  of  pre-averaging  K  and 
post-averaging  M.  In  particular,  numerous  ROCs  are  presented  in 
appendices  D,  E,  F,  and  G  for  the  cases  of  KM  =  4,  16,  64,  and 
256,  respectively.  The  particular  cases  run  are  indicated  in 
table  1;  a  total  of  124  ROCs  are  represented.  For  N  =  1,  since 
only  the  product  KM  matters,  only  one  ROC  need  be  presented  (the 
first  curve  in  each  appendix)  under  the  title  labeled  M  =  1. 


Table  1.  Cases  Run 


K 


M 


N 


4 

2 

1 


1 

2 

4 


1,2,4,8,16,32  Appendix  D 

1+15  ROCs 


16 

8 

4 

2 

1 


1 

2 

4 

8 

16 


1,2,4,8,16,32 


Appendix  E 
1+25  ROCs 


64 

32 

16 

8 

4 

2 

1 


1 

2 

4 

8 

16 

32 

64 


1,2,4,8,16,32 


Appendix  F 
1+35  ROCs 


256 

1 

128 

2 

64 

4 

32 

8 

16 

16 

8 

32 

4 

64 

2 

128 

1 

256 

1,2,4,8,16,32 


Appendix  G 
1+45  ROCs 
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NUMERICAL  APPROACH 


Whenever  the  CF  of  or-ing  output  v(t)  can  be  derived  in 
closed  form,  such  as  in  equations  (72),  (73),  (78),  or  (79),  it 
is  used  directly  in  equation  (35)  to  find  the  CF  of  decision 
variable  w(t).  However,  when  this  is  not  possible,  the  closed- 
form  results  for  the  or-ing  output  PDFs  given  in  equations  (32) 
and  (33)  are  used,  instead,  in  the  cascaded  FFT  procedure 
outlined  in  equations  (34)  and  (35).  In  the  former  case,  the 
fundamental  sampling  increment  is  A^,  starting  in  the  CF  domain, 
whereas  in  the  latter  case,  the  fundamental  sampling  increment  is 
A  starting  in  the  PDF  domain.  The  pertinent  numerical 
parameter  values  used  for  each  case  investigated  are  indicated  on 
the  appropriate  ROCs  in  appendices  D  through  G. 


REQUIRED  INPUT  SIGNAL-TO-NOISE  RATIOS 

The  numerous  ROCs  in  appendices  D  through  G  enable  a  user  to 

easily  investigate  and  determine  the  amount  of  system  input  SNR 

required  for  the  or-ing  device  to  achieve  a  specified  level  of 

detectability  performance,  in  terms  of  a  desired  false  alarm 

probability  Pf  and  detection  probability  Pd,  over  a  wide  range  of 

parameter  values,  in  an  effort  to  partially  condense  this 

voluminous  information,  a  standard  operating  point  (SOP)  is 

defined  here  as  Pf  =  IE-3,  Pd  =  0.5,  and  a  high-quality  operating 

point  (HOP)  is  defined  as  Pf  =  IE-6,  Pd  =  0.9.  The  corresponding 

2  2 

required  values  of  input  SNR  p(dB)  =  10  log(p)  =  10  1°9( ffs/ffn) j 


as  determined  from  the  ROCs  in  appendices  D  through  G,  are  listed 
in  tables  2  through  9.  (These  accurate  values  were  interpolated 
from  the  ROCs  while  the  false  alarm  and  detection  probability 
numbers,  and  P^,  were  still  resident  in  the  computer;  eyeball 
interpolation  from  the  plotted  ROCs  cannot  be  done  this 
accurately. ) 


These  input  SNR  requirements  are  plotted  in  figures  2  through 
9.  As  expected,  the  required  input  SNR  p  increases  roonotonically 
with  the  number  N  of  channels  or-ed  if  K  and  M  are  held  fixed. 
Also,  the  required  input  SNR  p  increases  monotonically  with  the 
amount  M  of  post-averaging  employed,  if  N  and  KM  are  held  fixed; 
alternatively,  the  required  input  SNR  p  decreases  with  the  amount 
K  of  pre-averaging  employed,  if  N  and  KM  are  held  fixed.  The 
exact  rates  can  be  determined  from  figures  2  through  9. 


Table  2.  Required  Input  SNR  p(dB)  for  P^  =  lE-3,  P^  =  0.5, 


KM  =  4, 

Gaussian 

Signal 

K 

M 

N  =  1 

N  =  2 

II 

Z 

CO 

11 

Z 

VO 

T— 1 

II 

Z 

N  =  32 

4 

1 

4.08 

4.46 

4.81 

5.13 

5.42 

5.70 

2 

2 

4.08 

4.66 

5.17 

5.62 

6.02 

6.38 

1 

4 

4.08 

4.82 

5.47 

6.05 

6.55 

7.00 

Table  3.  Required  Input 

SNR  p(dB) 

for  P^  = 

IE-6,  Pd  = 

0.9, 

KM  =  4, 

Gaussian 

Signal 

K 

M 

H 

II 

Z 

N  =  2 

ll 

55 

N  =  8 

N  =  16 

N  =  32 

4 

1 

10.50 

10.57 

10.84 

11.00 

11.16 

11.31 

2 

2 

10.50 

10.81 

11.10 

11.37 

11.59 

11.83 

1 

4 

10.50 

10.96 

11.37 

11.73 

12.08 

12.39 
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Table  4.  Required  Input  SNR  p(dB)  for  =  IE-3,  Pd  =  0.5, 

KM  =  16,  Gaussian  Signal 


K 

M 

N  =  1 

N  =  2 

N  =  4 

N  =  8 

N  =  16 

N  =  32 

16 

1 

-0.03 

0.31 

0.61 

0.89 

1.14 

1.38 

8 

2 

-0.03 

0.53 

1.00 

1.42 

1.79 

2.12 

4 

4 

-0.03 

0.73 

1.39 

1.95 

2.43 

2.86 

2 

8 

-0.03 

0.91 

1.72 

2.42 

3.02 

3.54 

1 

16 

-0.03 

1.03 

1.98 

2.81 

3.51 

4.12 

Table  5.  Required  Input 

SNR  p(dB) 

for  P^  = 

IE-6 ,  Pd  = 

0.9, 

KM  -  16 

,  Gaussian 

Signal 

K 

M 

tH 

li 

25 

N  =  2 

N  -  4 

CO 

II 

S 

N  =  16 

N  =  32 

16 

1 

4.51 

4.65 

4.79 

4.92 

5.05 

5.17 

8 

2 

4.51 

4.78 

5.03 

5.26 

5.47 

5.67 

4 

4 

4.51 

4.95 

5.34 

5.69 

6.01 

6.29 

2 

8 

4.51 

5.13 

5.67 

6.15 

6.57 

6.94 

1 

16 

4.51 

5.27 

5.95 

6.55 

7.08 

7.54 

Table  6.  Required  Input  SNR  p(dB)  for  Pf  =  IE-3,  Pd  =  0.5, 

KM  =  64,  Gaussian  Signal 


K 

M 

N  =  1 

N  =  2 

N  =  4 

GO 

II 

2 

N  =  16 

N  =  32 

64 

1 

-3.58 

-3.28 

-3.00 

-2.74 

-2.52 

-2.31 

32 

2 

-3.58 

-3.06 

-2.61 

-2.22 

-1.89 

-1.58 

16 

4 

-3.58 

-2.83 

-2.19 

-1.66 

-1.20 

-0.81 

8 

8 

-3.58 

-2.61 

-1.80 

-1.11 

-0.53 

-0.04 

4 

16 

-3.58 

-2.44 

-1.45 

-0.61 

0.09 

0.68 

2 

32 

-3.58 

-2.33 

-1.18 

-0.21 

0.61 

1.32 

1 

64 

-3.58 

-2.26 

-1.03 

0.07 

1.01 

1.83 

32 


Table  ’ 

7.  Required  Input 

KM  =  64 

SNR  p(dB) 

,  Gaussian 

for  P^  = 
Signal 

IE-6 ,  Pd  = 

o 

• 

VO 

K 

M 

rH 

II 

55 

N  =  2 

II 

00 

II 

Z 

N  =  16 

N  =  32 

64 

1 

0.12 

0.24 

0.36 

0.47 

0.57 

0.68 

32 

2 

0.12 

0.36 

0.57 

0.77 

0.96 

1.13 

16 

4 

0.12 

0.52 

0.88 

1.19 

1.48 

1.74 

8 

8 

0.12 

0.72 

1.24 

1.69 

2.08 

2.43 

4 

16 

0.12 

0.92 

1.60 

2.19 

2.70 

3.14 

2 

32 

0.12 

1.08 

1.91 

2.64 

3.27 

3.80 

1 

64 

0.12 

1.18 

2.15 

3.00 

3.73 

4.37 

Table  8.  Required  Input  SNR  p(dB)  for  Pf  =  IE-3,  Pd  =  0.5, 

KM  =  256,  Gaussian  Signal 


K 

M 

N  =  1 

N  =  2 

N  =  4 

N  =  8 

N  =  16 

N  -  32 

256 

1 

-6.87 

-6.57 

-6.32 

-6.07 

-5.87 

-5.65 

128 

2 

-6.87 

-6.37 

-5.94 

-5.57 

-5.24 

-4.96 

64 

4 

-6.87 

-6.12 

-5.51 

-5.00 

-4.56 

-4.18 

32 

8 

-6.87 

-5.89 

-5.08 

-4.42 

-3.87 

-3.39 

16 

16 

-6.87 

-5.67 

-4.69 

-3.86 

-3.18 

-2.60 

8 

32 

-6.87 

-5.52 

-4.35 

-3.36 

-2.55 

-1.86 

4 

64 

-6.87 

-5.42 

-4.10 

-2.97 

-2.01 

-1.21 

2 

128 

-6.87 

-5.36 

-3.96 

-2.69 

-1.60 

-0.68 

1 

256 

-6.87 

-5.36 

-3.90 

-2.53 

-1.33 

-0.30 

Table  9.  Required  Input  SNR  p(dB)  for  Pf  =  IE-6,  Pd  =  0.9, 

KM  =  256,  Gaussian  Signal 

K  M  N  =  1  N  =  2  N  =  4  N  =  8  N  =  16  N  =  32 

256  1  -3.57  -3.46  -3.36  -3.25  -3.15  -3.06 

128  2  -3.57  -3.36  -3.15  -2.98  -2.81  -2.64 

64  4  -3.57  -3.18  -2.86  -2.57  -2.31  -2.07 

32  8  -3.57  -2.99  -2.49  -2.07  -1.70  -1.39 

16  16  -3.57  -2.77  -2.10  -1.54  -1.06  -0.65 

8  32  -3.57  -2.57  -1.72  -1.01  -0.41  0.10 

4  64  -3.57  -2.42  -1.40  -0.54  0.19  0.80 

2  128  -3.57  -2.30  -1.16  -0.16  0.69  1.41 

1  256  -3.57  -2.25  -1.01  0.10  1.07  1.90 
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p(dB) 


M 


Figure  2.  Required  Input  SNR  for  Pf  =  IE  -  3,  Pd  =  0.5, 
EM  =  4,  Gaussian  Signal 


1  2  4 


M 

Figure  3.  Required  Input  SNR  far  Pf  =  IE  -  6,  Pd  = 
EM  =  4,  Gaussian  Signal 
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TRADEOFF  BETWEEN  PRE-  AND  POST-AVERAGING 


For  a  fixed  amount  of  total  time-bandwidth  product,  that  is, 
KM  constant,  it  is  possible  to  shift  some  of  the  pre-averaging  to 
post-averaging,  or  vice-versa.  For  example,  if  KM  =  16,  this 
could  be  accomplished  by  taking  K  =  16,  M  =  1,  which  corresponds 
to  no  post-averaging,  and  switching  to  K  =  1,  M  =  16,  which 
corresponds  to  no  pre-averaging.  Alternatively,  intermediate 
values  such  asK=4,  M=4,  could  be  used,  still  keeping 
KM  =  16. 


In  all  cases,  increasing  M  at  the  sake  of  K  requires  larger 
input  SNRs  to  maintain  the  same  operating  point;  see  figures  2 
through  9.  This  is  related  to  the  fact  that  the  optimum 
processor  corresponds  to  conducting  all  pre-averaging  before  any 
nonlinear  operations;  see  appendix  A.  For  example,  if  KM  =  16, 
in  order  to  maintain  the  SOP,  table  4  indicates  that  increases  of 
0.72  dB,  1.37  dB,  1.92  dB,  2.37  dB,  2.74  dB  are  required  for 
N  *  2,  4,  8,  16,  32,  respectively,  if  the  switch  is  made  from  all 
pre-averaging  to  all  post-averaging.  Alternatively,  for  the  HOP, 
the  corresponding  differences  in  table  5  are  smaller  (0.62  dB, 
1.16  dB,  1.63  dB,  2.03  dB,  2.37  dB  for  N  -  2,  4,  8,  16,  32, 
respectively) . 

For  a  HOP  like  P^  =  IE-6,  P^  =  0.9,  larger  input  SNRs  are 
naturally  required  to  achieve  this  level  of  performance.  At 
these  larger  input  SNRs,  the  or-ing  process  is  more  frequently 
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dominated  by  the  signal-bearing  channel.  Therefore,  as  the 
number  N  of  channels  or-ed  increases,  the  required  increases  in 
SNR  are  less  severe  with  increasing  N  than  for  a  lower-quality 
operating  point,  such  as  ■  IE-3,  =  0.5.  For  example,  for 

K  =  16,  M  =  1,  increasing  N  from  2  to  32  requires  a  1.07-dB  SNR 
increase  to  maintain  the  SOP,  whereas  only  an  additional  0.52  dB 
is  needed  at  the  HOP. 

A  related  observation  is  that,  for  a  given  configuration 
(fixed  K  and  M),  the  increase  in  SNR  required  to  maintain  the 
SOP,  as  the  amount  of  or-ing  increases,  is  more  severe  for  the 
larger  M  values.  Thus,  for  KM  =  16,  as  N  increases  from  2  to  32, 
a  1.07-dB  SNR  increase  (1.38  -  0.31)  suffices  to  maintain  the  SOP 
for  M  *  1,  whereas  a  3.09-dB  increase  (4.12  -  1.03)  is  required 
for  M  =  16.  The  corresponding  increases  at  the  HOP  are  0.52  dB 
(5.17  -  4.65)  for  M  =  1  versus  2.27  dB  (7.54  -  5.27)  for  M  =  16, 
as  N  increases  from  2  to  32.  The  behaviors  are  similar  for  the 
other  values  of  KM,  although  the  required  input  SNR  p  values  are 
smaller  as  KM  increases,  due  to  the  additional  observation  times 
allowed. 


39/(40  blank) 


EFFECT  OF  OVERLAPPING  DATA  AVERAGES 


Each  pre-averager  output  is  the  result  of  an  accumulation  of 
the  past  K  data  values  at  its  input.  Therefore,  each  adjacent 
pair  of  pre-averager  time  outputs  is  rather  highly  correlated, 
especially  for  large  K.  However,  the  or-ing  output  is  the  result 
of  a  comparison  of  N  independent  pre-averager  outputs,  with  only 
the  largest  surviving.  For  large  N,  this  competition  decreases 
the  correlation  of  successive  or-ing  time  outputs,  because  a 
different  channel  may  dominate  the  or-ing  decision  each  time 
sample  (see  appendix  H).  This  decrease  in  correlation  suggests 
that  the  or-ing  output  should  be  sampled  more  often  than  once 
every  K  samples,  which  corresponds  to  block  averaging  of 
non-overlapping  adjacent  blocks  of  data,  that  is,  skip  factor 
J  =  K. 


However,  there  are  two  conflicting  factors  that  tend  to  favor 
sampling  only  once  every  K  samples.  The  first  one  is  that,  for 
signal  present  with  a  fairly  reasonable  input  SNR,  the  same 
or-ing  channel  input  will  tend  to  dominate  the  decision,  thereby 
tending  to  keep  the  or-ing  output  correlation  at  larger  values, 
which  would  obviate  the  need  for  more  frequent  sampling. 

The  second  factor  is  that  when  the  sampling  interval  at  the 
or-ing  output  is  J  =  K,  the  effective  weighting  of  the  input  data 
is  uniform  over  the  total  observation  time  T.  This  uniform 
effective  weighting  would  be  optimum  for  detection  purposes  if 
there  were  no  or-ing,  N  =  1,  since  none  of  the  input  time  samples 
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are  preferred  over  any  other  time  samples  (see  appendix  A).  At 
the  other  extreme  of  sampling  every  successive  or-ing  time 
output,  the  effective  weighting  has  a  trapezoidal  shape,  peaking 
in  the  interior  of  the  observation  interval.  This  effective 
weighting  discounts  the  edge  samples  in  the  observation  interval 
in  terms  of  their  effect  on  the  decision  variable. 

Determining  which  effect  dominates  the  processor  performance 
can  be  ascertained  only  by  conducting  simulations  to  evaluate  the 
ROCs  for  a  variety  of  parameter  values,  including  the  amount  M  of 
post-averaging.  Each  of  the  three  results  in  tables  10,  11,  and 
12  (for  a  Gaussian  signal)  has  been  obtained  under  the 
simultaneous  conditions  of  keeping  the  pre-averaging  time  K 
fixed,  keeping  the  post-averaging  time  (M-l)J  fixed,  and, 
therefore,  keeping  the  total  observation  time  T  *  K  +  (M-l)J 
fixed.  The  three  cases  considered  correspond,  respectively,  to  a 
small  total  observation  time  T  =  24,  an  intermediate  value 
T  *  64,  and  a  large  value  T  =  256.  The  amount  of  or-ing  is 
varied  from  N  =  1  to  N  =  32. 

In  tables  10,  11,  and  12,  the  required  SNR  values  p(dB)  that 
are  not  enclosed  in  parentheses  correspond  to  the  SOP,  whereas 
the  values  in  parentheses  are  for  the  HOP.  The  latter  sets  are 
given  only  for  J  =  K  or  N  =  1,  both  of  which  special  cases  can  be 
evaluated  accurately  analytically.  The  remaining  SNR  values  in 
tables  10,  11,  and  12  were  obtained  by  simulations  using  one 
million  trials;  therefore,  estimation  of  the  required  input  SNR 
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to  achieve  P^  =  IE-6  is  not  possible  with  this  few  trials.  (The 
HOP  values  for  the  top  row,  N  =  1,  of  table  12  are  -3.57,  -3.53, 
-3.50,  -3.48,  -3.47  dB,  respectively.) 

Tables  10,  11,  and  12  are  plotted  in  figures  10,  11,  and  12. 
Perusal  of  these  results  reveals  that  the  losses  associated  with 
increased  or-ing,  N,  are  not  Significantly  ameliorated  by  using 
J  <  K,  when  the  amounts  of  pre-  and  post-averaging  are  kept 
constant.  For  example,  from  table  10,  with  a  Gaussian  signal, 

K  =  8,  N  =  32,  (M-l)J  =  16,  the  improvement  in  switching  from 
J  =  8  to  J  =  1  is  virtually  nil.  In  fact,  it  is  better  to  use 
J  =  4  for  this  particular  example  with  N  =  32,  that  is,  50% 
overlap,  for  which  the  improvement  is  approximately  0.1  dB  (1.44 
-  1.31)  relative  to  J  =  8.  The  reason  for  this  behavior  is  the 
effective  trapezoidal  weighting.  In  fact,  for  slight  or-ing, 
that  is,  small  N,  the  effective  trapezoidal  weighting  for  small  J 
actually  leads  to  losses  in  performance  relative  to  J  =  K. 

In  the  second  Gaussian-signal  example,  table  11,  with  K  =  4, 

N  =  32,  (M-l)J  =  60,  the  gain  at  the  SOP,  in  switching  from  J  =  K 
to  J  =  1,  is  only  0.3  dB  (0.68  -  0.36).  Smaller  N  values  achieve 
even  smaller  gains.  Finally,  the  long  observation  time  T  =  256, 
shown  in  table  12,  confirms  the  trends  anticipated  earlier.  In 
summary,  for  moderate  amounts  of  or-ing,  overlapped  averaging  is 
not  worthwhile.  On  the  other  hand,  for  a  significant  amount  of 
or-ing,  N  >>  1,  the  gains  would  be  more  substantial,  and  the  use 
of  J  <  K,  say  50%  overlap,  might  be  worth  re-considering. 
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Table  10.  Required  Input  SNR  p(dB)  for  =  IE-3,  =  0.5, 


K  = 

8,  (M-l)J 

=  16, 

II 

to 

* 

Gaussian  Signal 

N 

M=3 

CO 

il 

M=5 , 

J=4 

M=9,  J=2 

M=17 

1 

-1.10 

(3.12) 

-0.86 

(3.43) 

-0.68  (3.66) 

-0.56 

2 

-0.43 

(3.47) 

-0.32 

-0.18 

-0.08 

4 

0.14 

(3.79) 

0.14 

0.25 

0.33 

8 

0.63 

(4.08) 

0.57 

0.64 

0.71 

16 

1.06 

(4.34) 

0.97 

1.01 

1.07 

32 

1.44 

(4.58) 

1.31 

1.34 

1.38 

Table  11.  Required  Input  SNR  p(dB)  for  P^  =  IE-3,  P^  =  0.5, 


K  = 

=  4,  (M-l)J  = 

=  60, 

T  = 

64, 

daussian  Signal 

N 

M=16 , 

J= 

=4 

M 

=31 , 

,  J=2 

M=61 , 

J=1 

1 

-3.58 

(0, 

.12) 

-3 

.54 

(0.17) 

-3.51 

(0.20) 

2 

-2.44 

(0, 

.92) 

-2 

.50 

-2.51 

4 

-1.45 

(1. 

.60) 

-1 

.60 

-1.61 

8 

-0.61 

(2. 

.19) 

-0 

.82 

-0.86 

16 

0.09 

(2. 

.70) 

-0 

.13 

-0.19 

32 

0.68 

(3. 

.14) 

0 

.41 

0.36 

Table  12.  Required  Input  SNR  p(dB)  for  Pf  =  IE-3,  P^  =  0.5, 
K  =  16,  (M-l)J  =  240,  T  =  256,  Gaussian  Signal 


N 

M=16 

,  J=1 6 

M=31 , J=8 

M-61, J=4 

M=121 , J=2 

M=241 , J=1 

1 

-6.87 

(-3.57) 

-6.83 

-6.80 

-6.79 

-6.78 

2 

-5.67 

(-2.77) 

-5.77 

-5.78 

-5.77 

-5.77 

4 

-4.69 

(-2.10) 

-4.89 

-4.92 

-4.91 

-4.92 

8 

-3.86 

(-1.54) 

-4.10 

-4.16 

-4.18 

-4.17 

16 

-3.18 

(-1.06) 

-3.47 

-3.53 

-3.55 

-3.55 

32 

-2.60 

(-0.65) 

-2.93 

-3.01 

-3.03 

-3.04 
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Figure  12.  Required  Input  SNR  for  Pf  =  IE  -  3,  Pd  = 
K  =  16,  (M  -  1)J  =  240,  T  =  256,  Gaussian  Signal 


SUMMARY 


For  three  different  signal  models  in  Gaussian  noise,  closed 
forms  have  been  derived  for  the  PDFs  at  the  output  of  a  pre¬ 
averager  followed  by  an  or-ing  operation.  In  the  case  of  a 
random  Gaussian  signal,  these  forms  have  been  numerically 
transformed  to  yield  accurate  results  for  the  EDFs  at  the  output 
of  a  post-averager,  under  both  hypotheses  Hq  and  H^,  for  numerous 
values  of  the  parameters  of  the  complete  processor.  These 
results  enable  investigation  of  the  or-ing  processor  for  false 
alarm  probabilities  in  the  range  of  IE-6  and  smaller;  there  is 
no  need  to  resort  to  lengthy  simulations. 

Numerous  ROCs  have  been  generated,  which  enable  a  user  or 
system  designer  to  quickly  assess  the  losses  to  be  expected  from 
employing  or-ing  in  a  processor.  Also,  quantitative  evaluation 
of  the  tradeoffs  between  pre-averaging  versus  post-averaging  has 
been  conducted.  Finally,  since  the  enclosed  tabulations  will 
undoubtedly  not  cover  all  cases  of  practical  interest,  a  MATLAB 
program  for  evaluation  of  the  ROCs  is  presented  in  appendix  I  for 
additional  evaluation  and  tabulation. 

The  possibility  of  employing  overlapping  data  in  the  post¬ 
averager  has  been  investigated  by  simulation  and  found  to  be 
insignificant,  at  least  for  moderate  values  of  N,  the  number  of 
channels  or-ed.  For  very  large  values  of  N,  this  conclusion  may 
require  reconsideration. 
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Numerical  results  have  been  presented  for  KM  =  4,  16,  64,  and 
256.  For  completeness,  the  special  case  of  KM  =  1  is  presented 
here.  The  is  given  by  a  combination  of  equations  (37)  and 
(68),  while  is  given  by  equations  (36)  and  (69),  with 
K  =  1.  For  given  values  of  P^  and  N,  the  required  threshold 
value  u  can  be  solved.  Then,  for  a  given  P^  and  the  above 
solution  for  u,  the  required  input  SNR  p  can  be  found,  using 
a  =  1/(1  +  p),  in  the  explicit  form  for  a  Gaussian  signal  as 


_  ln(  1  -  <x)  _ 

p  "  ln( 1  -  ap) 


1/N 


(83) 


This  relation  is  plotted  in  figure  13  as  a  function  of  the  amount 
of  or-ing  N  for  both  the  SOP  and  the  HOP. 


The  two  alternative  cases  of  a  phase-incoherent  signal  and  a 
deterministic  signal  in  additive  Gaussian  noise  are  currently 
under  investigation.  Tabulation  of  the  corresponding  ROCs  and 
tradeoffs  between  pre-  and  post-averaging  will  be  conducted  in  a 
similar  manner  and  for  the  same  parameter  values  as  accomplished 
here . 


N 


Figure  13.  Required  Input  SNR  for  K  =  1,  M  =  1 
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APPENDIX  A  -  OPTIMUM  PROCESSING  OP  TOTAL  DATA  SET 


The  observed  data  in  the  n-th  channel  are  xn(t)  for  1  <  t  <  T. 
There  are  N  channels  giving  rise  to  the  total  data  set 

{x  (t)}  s  x  ( t )  for  1  <  n  <  N  ,  1  <  t  <  T  .  ( A-l ) 

n  n 

All  N  T  of  these  samples  are  Statistically  independent  of  each 
other.  Some  of  the  following  analysis  was  presented  earlier  in 
reference  3,  appendix  E. 

It  is  desired  to  determine  the  a  posteriori  probability  (APP) 
that  channel  number  j  contains  a  signal,  namely, 

Pr(j)  p.({x  (t)}) 

Pr(j  |{xn(t)>)  =  - p  -  for  0  <  j  <  N  .  (A-2) 

The  case  j  =  0  corresponds  to  no  signal  present.  For  j  >  1, 

Pr(j)  is  the  a  priori  probability  of  the  signal  occupying  channel 
j  (during  the  entire  observation  interval  T),  and  p j ( * )  is  the 
PDF  of  total  observation  (x  (t)}  when  the  signal  is  in  channel  j. 
Denominator  Pc(*)  is  the  combined  PDF  of  the  total  data,  namely, 
the  sum  of  the  numerator  terms  from  j  =  0  to  N. 

Let  p  denote  the  PDF  of  individual  sample  xn(t)  for  signal 
absent,  and  let  £  denote  the  PDF  of  xn(t)  for  signal  present. 
Then,  the  j  =  0  PDF  in  the  numerator  of  equation  (A-2)  is 

N  T 

p  ({x  (t)})  =  TT  TT  P(x  ( t )  )  ,  ( A-3 ) 

u  n  n=l  t=l  n 


A-l 


while  the  other  terms  are 


Pj<{*n(t)})  = 


r  T  'j  N  T 

TT  E(x.(t) )  1  TT  "TT 
lt-1  3  J  n=l  t-1 


p ( Xn ( t ) ) 


n*j 


-  P0({xn(t)}) 


T 

t=l 


"~r  l(  x .  ( t ) ) 


for  1  <  j  <  N  , 


(A-4) 


where  auxiliary  function 


L(u) 


E(u) 

PUT  * 


( A-5 ) 


The  APPs  in  equation  (A-2)  are 

Pr ( 0 | • )  =  Pr ( 0 )  P0 ( * )/Pc ( ' )  ' 

T 

Pr  ( j  I  • )  =  Pr  ( j  )  pn(-)/p„(*)  TT  Mx.(t))  for  1  <  j  <  N  .  (A-6) 

t=l  J 

Optimum  processing  corresponds  to  selection  of  the  channel  with 
the  largest  APP.  Upon  canceling  common  factors,  this  rule  leads 
to  consideration  of  the  maximum  of  N+l  quantities  according  to 

(  T  T  'i 

max  jpr  ( 0 ) ,  Pr(l)  TT  Mx^t)),  /  Pr(N)  TT  L(xN(t>)|  •  <A"7) 

If  the  zeroth-order  term  is  largest,  a  declaration  is  made  that 
no  signal  is  present.  If  the  j-th  term  is  largest,  j  >  0,  a 
signal  is  declared  present  in  channel  j. 

An  alternative  processing  procedure  is  obtained  when  channel 
identification  is  not  of  interest;  rather,  only  the  presence  or 


A-2 


absence  of  a  signal  is  relevant.  Then,  for  equal  a  priori 
probabilities,  Pr(j)  =  Pr(l)  for  1  <  j  <  N,  the  APP  of  signal 
present  somewhere  is,  upon  use  of  equations  (A-2)  and  (A-4), 


Pr(l) 

Pc<*> 


N 

u 

j=i 


Pj(  •) 


Pr(l) 


Po<->  3, 

n  f  L— 4 


T 

FT  L(x.(t))  . 

t=l  3 


(A-8) 


The  alternative  APP  of  signal  absent  is 


Pr(0) 


P0(‘) 

Pc(,)  * 


( A-9 ) 


Picking  the  largest  of  these  two  quantities  is  tantamount  to  a 
threshold  comparison  according  to 


NT 

ZU  FT  Mx  (t))  <  v  .  ( A-10 ) 

n=l  t=l  n  < 


The  two  alternative  test  statistics  (A-7)  and  (A-10)  will  now  be 
calculated  for  the  three  signal  models  of  interest  here. 


COHERENT  (DETERMINISTIC)  SIGNAL 

p(u)  =  (2  Jiff2)-5*  exp(-  ,  £(u)  =  (  2ji<j2  )-J*  exp  j  , 

f(S)  «  exp( i £m  -  £2a2/2)  , 

L(u)  -  exp(2|  -  |  j  ,  d  =  ~  .  (A— 11 ) 

O 

The  CF  corresponding  to  PDF  £(u)  is  f  ( U  .  For  1  <  n  <  N,  the 


n-th  data  processing  term  in  equation  (A-7)  is 


T  T 

]  [  L(x  (t))  =  exp(-Td2/2)  exp  YU  xn(t))  *  (A-12) 

t=l  n  ^  t=l  n  * 

Therefore,  the  optimum  test  statistic  takes  the  form  of  threshold 
comparison 


T 

■  E 

1  <n<N  lt=l 


max 


I: 


xn(t> 


> 

< 


(A-13 ) 


Physically,  each  channel  output  is  summed  over  the  entire 
observation  interval  T,  and  the  largest  channel  sum  is  compared 
with  a  threshold  for  declaration  of  signal  presence  and  channel 
number.  If  none  exceed  the  threshold,  signal  is  declared  absent. 


On  the  other  hand,  when  channel  identity  is  irrelevant, 
equations  (A-10)  and  (A-12)  yield  test 

EZ  exp(f  C  x  (t))  >  v  .  (A-14) 

Again,  each  channel  is  summed  in  time,  but  must  now  be  scaled  by 
d/cr,  exponentiated,  and  summed  over  all  channels  prior  to 
threshold  comparison.  For  reasonable  input  SNRs  d,  the  outer  sum 
on  n  will  be  dominated  by  the  exponential  with  the  largest  inner 
sum  on  t;  in  that  case,  test  (A-14)  substantially  reduces  to  test 
(A-13).  This  is  fortunate,  because  scaling  d/a  will  not  likely 
be  known  in  practice.  Thus,  test  (A-13)  is  essentially  the 
optimum  processor  for  both  cases  of  whether  channel 
identification  is  desired  or  not. 
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PHASE-INCOHERENT  SIGNAL 


p(u)  =  ij  Mp(-f)  -  E(»)  -  ^2  ®*p(-t  -p)  Io(2("U>'1)  ' 

<T  O’  O  C 

2 

f(U  =  - - - J  exp fp  U<T  j)  ,  xTTtT  =  a2 ( l  +  p)  , 

1  -  Ua  k  1  -  Her  ' 


h  2 

L (u )  *  exp(-p)  I0(2(^u)  )  ;  P  =  ,  u  >  0  .  (A-15) 


For  1  <  n  <  N,  the  n-th  data  processing  term  in  equation  (A-7)  is 
T  t 

FT  L(xn(t))  =  exp(-Tp)  77  1 {xn(  t)  J*5)  .  ( A- 1 6  ) 


It  is  important  to  notice  that  the  arguments  of  the  Bessel 

function  IQ  have  not  undergone  any  summations,  either  on  t  or  n. 

(This  is  in  contrast  to  test  (A-14),  where  the  arguments  of  exp 

are  summed  on  t  before  exponentiation. )  Therefore,  these 

arguments  in  equation  (A-16)  are  typically  small,  allowing  the 

2 

approximation  lQ(z)  =  exp(z  /4)  for  small  z,  since 


2  ,  2..  .2 

IQ(z)  =  1  +  f-  +  +  •••  , 

o  2  ,  2..  x2 

exp( z2/4 )  =  1  +  f-  +  (z  ^4)  +  •••  . 


(A-17) 


This  leads  to  an  approximation  for  equation  (A-16),  namely, 


T  T 

exp(-Tp)  TT  exp(£2  xn(t0  "  exP("Tp)  exp^XZl  xn(t))  •  (A”18 
t=l  a  va  t=l 


Therefore,  the  optimum  test  statistic  takes  the  form  of 
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threshold  comparison 


( A-19 ) 


Although  this  approximate  test  has  a  form  identical  to  equation 
(A-13),  xR(t)  here  is  a  squared  envelope,  whereas  it  was  an 
amplitude  quantity  previously;  see  equation  (A-15)  versus 
equation  (A-ll). 


On  the  other  hand,  processor  (A-10),  with  the  help  of 
equation  (A-18),  yields  the  alternative  test 

LZ  «*p(£5  C  V1’)  <  v  .  I*"2®) 

n=l  a  t=l 

2 

The  N  temporal  sums  are  scaled  by  p/a  ,  exponentiated,  and  then 
summed  over  all  channels.  For  reasonable  input  SNRs  p,  the  n-sum 
is  dominated  by  the  exponential  with  the  largest  t-sum;  thus, 
this  test  is  essentially  equivalent  to  test  (A-19).  Therefore, 
test  (A-19)  is  virtually  the  optimum  processor  under  both  cases 
of  channel  identification. 


ACCESSIBILITY  TO  IN-PHASE  AND  QUADRATURE  COMPONENTS 

The  statistical  model  in  equation  (A-15)  corresponds  to 
sampling  the  squared  envelope  of  the  received  data.  If,  instead, 
the  in-phase  and  quadrature  components  are  accessible,  input  data 
value  *n(t)  in  the  signal  channel  should  be  replaced  by  the  two 
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RVS 


ln ( t )  =  A(t)  cosO(t)  +  g^t)  , 

Qn(t)  =  A(t)  sin0( t )  +  gg(t)  ,  (A-21) 

where  amplitude  A(t)  is  non-random,  RV  0(t)  is  uniformly 
distributed  over  2n,  and  g^(t)  and  gq(t)  are  independent, 
identically  distributed  Gaussian  RVs  with  zero  mean  and  variance 
<y2.  The  joint  OF  of  RVs  In(t)  and  Qn(t)  is 

2 

fU,n)  =  exp[iUn(t)+inQn(t)]  =  J0(A(t)U2+h2)V)  exp(-  f-U2+h2)). 

( A-22 ) 


From  this  point  on,  it  is  presumed  that  amplitude  A(t)  is 
constant,  that  is,  A(t)  =  A  for  all  t  in  the  observation 
interval.  The  joint  PDF  corresponding  to  equation  (A-22)  is  then 


£(I,Q) 


(A-23) 


Since  this  expression  doesn't  factor  in  I  and  Q,  the  in-phase  and 

quadrature  components  are  statistically  dependent  on,  each  other. 

2  2 

However,  equation  (A-23)  depends  only  on  the  combination  I  +  Q  , 
and  not  on  I  or  Q  separately;  any  optimum  processor  derived  using 
the  joint  PDF  (A-23)  must  process  the  accessible  data  In(t)  and 
Qn(t)  in  the  form  I2(t)  +  Q2(t).  In  fact,  if  one  identifies 
original  data  *n(t)  =  +  Q2(t)]/2,  the  resulting  PDF  p(u) 

for  *n(t),  as  obtained  from  equation  (A-23),  is  precisely 
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b  2 

equation  (A-15),  where  p  =  SjA  /a  .  The  end  result  is  that 
accessibility  to  in-phase  and  quadrature  data  values  (A-21)  would 
yield  the  same  optimum  processor  as  that  required  to  use  squared 
envelopes  {xn(t)}  directly. 


RANDOM  GAUSSIAN  SIGNAL 


p(u)  -  12  exp  ^  2)  '  2<u)  -  2  2  exp(  2  2)  ' 


a  +  a 


a  +  a 


f(0  = 


1  -  U(t2(1  +  p) 


*n(t>  «  » (1  +  p)  r 


L(u)  -  exP  ?  P  “  iW  ,  u  >  0  . 


( A-24 ) 


The  n-th  data  processing  term  for  1  <  n  <  N  in  equation  (A-7)  is 


t=l 


L ( xn ( t ) ) 


=  (1  +  p) 


-T 


exp 


(r 


+  P 


T 

t=i 


‘n(t0 


( A-25 ) 


Therefore,  the  optimum  test  statistic  takes  the  form  of  threshold 
comparison 


T 

■  E 

l<n<N  U=1 


max 


I: 


xn(t> 


}<" 


( A-26 ) 


Although  this  test  has  a  form  identical  to  test  (A-13),  xn(t)  is 
an  envelope-squared  quantity  here,  whereas  it  was  an  amplitude 
quantity  previously;  see  equation  (A-24)  versus  equation  (A-ll). 
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On  the  other  hand,  processor  (A-10)  yields  the  alternative 


test 

N  T 

1 1]  exp  fc-P— -  ^2  xn(t0  <  V  *  (A-27) 

n=l  u  +  p  <r  t-1  n  J  < 

The  N  temporal  sums  are  scaled,  exponentiated,  and  then  summed 
over  all  channels.  For  reasonable  input  SNRs  p,  the  n-sum  is 
dominated  by  the  exponential  with  the  largest  t-sum,  and  this 
test  is  essentially  equivalent  to  test  (A-26).  Thus,  test  (A-26) 
is  virtually  the  optimum  processor  under  both  cases  of  channel 
identification. 


ACCESSIBILITY  TO  IN-PHASE  AND  QUADRATURE  COMPONENTS 

The  statistical  model  in  equation  (A-24)  corresponds  to 
sampling  the  squared  envelope  of  the  received  data.  If,  instead, 
the  in-phase  and  quadrature  components  are  accessible,  input  data 
value  *n(t)  in  the  signal  channel  should  be  replaced  by  the  two 
RVs 

In(t)  =  ai(t)  +  9iU)  ' 

Qn<t)  =  aq(t)  +  9q(t)  '  (A-28) 

where  all  the  quantities  are  independent  zero-mean  Gaussian  RVs, 

2 

with  signal  samples  and  ag  having  variance  c  ,  while  noise 

2 

samples  and  have  variance  a  . 


The  joint  PDF  for  the  RVs  in  equation  (A-28)  is  obviously 


£(I,Q) 


1 

2n(  <r2+a2 ) 


I2  ±  Q2) 

2  ( <r2+ff2 )  ' 


(A-29) 


Thus,  In(t)  and  Qn(t)  in  equation  (A-28)  are  independent  RVs. 

More  importantly,  PDF  (A-29)  depends  only  on  the  combination 
2  2 

I  +  Q  ,  and  not  on  I  or  Q  separately;  any  optimum  processor 
derived  using  joint  PDF  (A-29)  must  process  the  accessible  data 
In(t)  and  Qn(t)  in  the  form  l£(t)  +  Q2(t).  In  fact,  if  one 
identifies  original  data  value  xn(t)  =  tl2(t)  +  Q^(t)]/2,  the 

resulting  PDF  jd(u)  for  xn(t),  as  obtained  from  equation  (A-29), 

2  2 

is  precisely  equation  (A-24),  where  p  =  a  /a  .  The  end  result  is 
that  accessibility  to  in-phase  and  quadrature  data  values  (A-28) 
would  yield  the  same  optimum  processor  as  that  required  to  use 
squared  envelopes  (xn(t)}  directly. 
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APPENDIX  B  -  ON  THE  USE  OF  CASCADED  FAST 
FOURIER  TRANSFORMS  FOR  DISTRIBUTION  CALCULATION 

The  PDF  of  or-ing  output  v(t)  is  often  available  in  closed 
form,  albeit  complicated,  for  signal  present  as  well  as  signal 
absent;  see  equations  (32)  and  (33).  For  independent  inputs  to 
the  post-averager,  that  is,  skip  factor  J  =  K,  the  CF  of  the 
post-averager  output  is  the  M-th  power  of  its  input  CF.  In  order 
to  find  the  post-averager  output  PDF  and  EDF,  its  input  PDF  must 
be  Fourier  transformed  to  find  the  corresponding  input  CF,  then 
raised  to  the  M-th  power,  followed  by  another  Fourier  transform 
back  to  the  density  or  distribution  domains.  This  cascaded  FFT 
operation  requires  care  in  its  evaluation  in  order  to  control 
aliasing  errors. 

Let  p(u)  be  the  post-averager  input  PDF  of  interest.  If  this 
function  is  sampled  at  increment  Au,  and  Fourier  transformed,  the 
result  is  an  approximation  to  desired  CF 

f(S)  =  J  du  exp(i^u)  p(u)  ,  (B-l) 

namely,  continuous  function 

fU)  =  Au  2  exp(UnAu)  p(nAu)  =  £  f  f  l  -  m  §*)  .  (B-2) 

n  m  u 

This  aliased  CF  f(£)  has  period  2ji/Au  in  £.  In  order  to  minimize 
the  aliasing  inherent  in  f(£),  sampling  increment  Au  must  be 
taken  small  enough  to  separate  the  lobes  in  equation  (B-2)  and 
avoid  significant  overlap. 


B-l 


Now,  the  aliased  CF  I(£)  will  be  evaluated  at  special 
selected  intervals,  namely, 

f  “  Au  E  exp(i2nkn/N1)  p(nAu)  for  0  <  k  <  Nj-1  ,  (B-3) 

which  covers  one  period  of  2(5).  and  k  are  integers. 

The  samples  of  p(u)  must  be  taken  far  enough  out  in  u  that 
the  two  neglected  tails  yield  insignificant  truncation  errors. 

If  the  number  of  u  samples  exceeds  ,  the  sample  at  u  =  nAu  can 
simply  be  prealiased  into  bin  (n  MODULO  N^),  since  the  factor 
exp( i2nkn/N1 )  treats  all  p(u)  samples  separated  by  equally. 

Thus,  the  truncation  error  in  u  can  be  controlled  to  any  desired 
degree  by  taking  enough  samples  of  PDF  p(u)  on  its  tails.  Equa¬ 
tion  (B-3)  can  be  accomplished  numerically  by  an  N^point  FFT. 

Proper  selection  of  sampling  increment  Au  can  be  achieved  by 
looking  at  a  plot  of  the  magnitude  of  equation  (B-3);  see  figure 
B— 1 .  The  right  skirt  of  the  desired  lobe  centered  at  5  =  0 


Figure  B-l.  Magnitude  of  Approximate  CF  f 
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(k  *  0)  will  cross  the  left  skirt  of  the  aliasing  lobe  centered 
at  £=2ji/Au  (k  =  Nx)  at  £  =  n/Au  (k  =  N1/2).  However,  it  is  not 
necessary  that  |f(£)|  be  extremely  small  at  this  crossover  at 
k  =  N^/2 .  Rather,  it  is  necessary  that  the  extrapolation  of  the 
left  skirt  of  the  aliasing  lobe  be  very  small  at  £  =  0  (k  =  0); 
this  result  can  be  visually  observed.  If  origin  value  V  of  the 
extrapolation  (dashed  curve)  is  not  sufficiently  small,  then  Au 
must  be  decreased  and  the  procedure  repeated  before  going  on. 


The  magnitude  of  the  absolute  error  of  equation  (B-2), 
namely,  |2(U  -  f(£)|  for  0  <  Z,  <  n/Au,  is  typically  (but  not 
always)  less  than  the  error  at  £  =  0.  Thus,  estimating  and 
controlling  this  origin  error  (by  means  of  Ay)  is  an  essential 
step.  This  error  can  be  measured  quantitatively  by  printing  out 
2(0).  Its  discrepancy  from  the  desired  value  f(0)  =  1  can  serve 
as  an  indicator  whether  to  increase  or  decrease  Au  before 
proceeding  any  further. 


Notice  that  FFT  size  has  nothing  to  do  with  the  accuracy 
of  approximation  2  in  equation  (B-3).  Rather,  merely  sets  the 
increment 


2n 


(B-4) 


at  which  samples  of  2($)  are  taken  in  the  £  domain  according  to 
equation  (B-3).  However,  N1  will  become  a  very  relevant 
parameter  in  the  remaining  part  of  the  cascade  FFT  procedure. 
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Once  a  satisfactory  value  of  Au  has  been  ascertained  (by 
trial  and  error)  in  order  to  realize  a  sufficiently  small  error 
in  f (0 ) ,  the  approximation  to  the  post-averager  output  CF  is 
computed  at  the  same  selected  intervals  according  to 

*(§^)  *  £or  0  * k  s  “i-1  •  (B-5) 

(The  values  for  k  >  N^/2  correspond  to  negative  $  values.) 
Although  the  relative  error  of  £  in  the  neighborhood  of  k  =  N^/2 
is  100%,  these  values  are  suppressed  to  insignificance  by  the 
M-th  power.  However,  the  small  error  in  f  at  k  =  0  is  magnified 
by  a  factor  of  approximately  M,  and  shows  up  in  approximation  f. 
Therefore,  large  values  of  M  may  require  that  Au  be  decreased 
even  further  in  anticipation  of  this  inherent  error  magnification 
caused  by  the  M-th  power  in  equation  (B-5). 

The  approximation  to  the  post-averager  output  PDF  is  given  by 
continuous  function  (inverse  Fourier  transform) 

N,-l 

p(u)  S  3^-  A-  YU  exp(-iukA-)  f(kA-)  .  (B-6) 

ZJl  *  k=0  **  ** 

This  function  has  period  2n/A^  =  N^Au  *n  u*  order  that  the 

inherent  aliasing  lobes  in  jd(u)  not  severely  contaminate  the 
approximation,  period  N^AU  must  be  made  large  enough  that  the 
aliasing  lobes  of  2<u)  do  not  overlap  significantly.  This  can  be 
ascertained  by  computing  approximation  £(u)  at  the  selected 
points  Uj  =  j  $u  h  j  2n/(N2A^),  getting 
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for  0  <  j  <  N2-l,(B-7) 


which  covers  one  period  of  p(u).  This  can  be  done  by  an  N2~point 
FFT.  If  N2  <  N^,  simply  prealias  f ( 2nk/(N1Au) )  into  bin 
(k  MODULO  N2),  since  the  exponential  factor  exp(-i2njk/N2 )  treats 
all  f  samples  separated  by  N2  equally.  The  u-increment  in 
equation  (B-7)  is 


5u  - 


2n 

N2A$ 


(B-8) 


in  terms  of  the  original  sampling  increment  Ay  on  input  PDF  p(u). 
Integer  N2  can  be  taken  larger  or  smaller  than  Nj,  since  it 
merely  controls  the  increment  at  which  approximation  p  is 
computed  in  equation  (B-7). 


In  the  special  case  N2  =  ,  then  5y  =  Au,  which  is  the  same 

sampling  increment  used  on  input  PDF  p(u).  Since  the  post¬ 
averager  output  is  the  sum  of  M  variables,  the  standard  deviation 
is  M*5  larger  than  at  the  summer  input,  making  the  output  PDF 
spread  over  a  larger  range.  If  Au  was  small  enough  to  track 
p(u),  then  Su  =  Ay  is  probably  fine  enough  to  track  p(u). 

Integer  N2  has  no  effect  on  the  accuracy  of  approximation  p. 
However,  integer  N1  in  FFT  (B-3)  now  has  a  significant  effect  on 
the  accuracy  of  the  calculated  PDF  at  the  post-averager  output, 
namely,  equation  (B-6).  Period  N-j^  of  p(u)  must  be  larger  than 
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the  extent  of  the  actual  unknown  output  PDF.  This  aliasing 
effect  in  jo(u)  can  be  observed  by  plotting  equation  (B-7),  which 
covers  one  period,  and  looking  for  overlap  of  the  PDF  tails 
before  they  have  sufficiently  decayed.  Since  the  period  of  jo(u) 
is  N^^u/  it  can  be  increased  only  by  increasing  ,  because 
has  already  been  selected  as  large  as  possible  without  causing 
significant  aliasing  in  f(£). 

A  limitation  could  arise  at  this  point.  If  Au  had  to  be 
chosen  very  small  because  input  CF  f(£)  was  a  very  slowly 
decaying  CF  (see  equation  (B-2)),  may  have  to  be  so  large  that 
sufficient  storage  is  not  available  to  conduct  N^-point  FFT 
(B-3).  (For  example,  PDFs  (56)  and  (71)  behave  as  u^  ^  as 
u  ->  0+ . )  One  possibility  is  then  to  break  p(u)  into  the  sum  of 
two  parts,  the  first  being  a  simple  analytic  function  that  takes 
on  the  highest  singularities  of  p(u)  and  has  a  known  closed-form 
Fourier  transform.  The  second  part,  which  is  smoother  in  u,  will 
have  a  Fourier  transform  that  decays  faster  and  can  profitably  be 
numerically  transformed  via  an  FFT.  The  sum  of  the  two  Fourier 
transforms  is  the  approximation  to  the  desired  CF  f(£). 

For  present  purposes,  where  false  alarm  and  detection 
probabilities  are  of  interest,  one  does  not  actually  compute  PDF 
£  in  equations  (B-6)  or  (B-7),  but  instead  proceeds  directly  from 
CF  f  to  the  corresponding  EDF.  This  accurate  efficient  procedure 
is  given  in  reference  9.  Nevertheless,  the  criterion  above  for 
choosing  FFT  size  N^  is  still  directly  relevant. 


B-6 


APPENDIX  C  -  MOMENTS  OF  CHARACTERISTIC  FUNCTION  (72) 


The  CF  of  or-ing  RV  v(t)  for  a  random  Gaussian  signal  and 
K  =  1  is 


£vi'E'  -  (1>N-1  +  (i-u>NJ  - 


a  =  1/(1  +  p)  ,  (C-l) 


under  H^;  see  equation  (72).  In  order  to  determine  the  moments 
of  v(t),  consider  the  following  development: 


(b-U),  =  fT(b-i^j-l)  , 
J  j=l 


(C-2) 


giving 


ln[ (b-i^)j] 


C  ln(b+j-l-i£) 

j-1 


g  ln(b+ j-l )  +  g  i«(i  -  sifn) 


J  ® 


ln[(b)a] 


(C-3) 


where 


3k( J,b) 


J 

-  n 


for  k  >  1 


j=l  (b+j-1) 


(C-4) 


From  equation  (C-3),  there  follows 


(b-iOj  -  (bjj  exp C  Bku,b>) 


( C-5 ) 
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Now,  identifying  J  -*  N,  b  ->  a  in  equation  (C-5),  one  has 


(1) 


i£ 


■  (%  eXP(5  <VN'a0  •  (c-6) 


N-l  (a-i£) 


N 


On  the  other  hand,  identifying  J  ->  N-l,  b  -»  1  in  equation  (C-5), 


( C-7 ) 


Now,  CF  in  equation  (C-l)  is  given  by  the  sum  of  equations 

( C-6 )  and  (C-7). 

To  simplify  the  two  relations  in  equations  (C-6)  and  (C-7), 
the  following  expansion  (reference  13,  appendix  A)  is  used: 


CD  ]{  CD 

exp  (J  ]  (l^  Pk(J,b)j  =  YZ  (i^)k  Yk(J,b)  , 

k— 1  k  ~  0 


( C-8 ) 


where  coefficients 


i  k 

Y0(J,b)  =  1  ,  yk(J,b)  =  £  C  ( J,b)  Yk_j(J,b)  for  k  >  1.  (C-9) 

Use  of  equation  (C-8)  on  equations  (C-6)  and  (C-7)  gives  CF 
(C-l)  as 


(  1  )  ®  CD 

fvim  =  — T~  ZZ  ( i^)k+1  rk(N,a)  +  JZ  <U)k  Yk(N-l,l).  (C-10) 


(a) 


N  k=0 


k=0 


The  moments  { /#  ( k ) }  of  or-ing  output  v(t)  are  then 
immediately  available  as 
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The  cumulants,  X^k),  of  v(t)  can  be  readily  found  by  using  the 
recursion  in  reference  13,  equation  ( C— 7 ) s 


vk>  -  "v(k)  -  C  (V)  Vk-j)  "v(j)  £or  k  >  1  •  (C 

Defining  normalized  statistics 

£v(k)  =  //v( k)/kl  ,  x^k)  =  x^kJ/U-l)!  ,  (C 

equation  (C-12)  becomes  (also,  see  reference  13,  equation  (C- 

k-1 

X^k)  =  k  £v(k)  -  YZ  Vk"j)  £v(j)  for  k  >  1  .  (C 

In  summary,  if  moments  {//v(k)}  and  cumulants  (X^k)}  of 
or-ing  output  v(t),  with  signal  present  and  K  =  1,  for  order 
k  =  1  to  J  are  of  interest,  the  following  calculations  are 
required,  with  a  =  1/(1  +  p): 

from  equation  (C-4),  evaluate  {^(I^a)}  up  to  k  =  J-l; 

from  equation  (C-4),  evaluate  {{^(11-1,1)}  up  to  k  =  J ; 

from  equation  (C-9),  evaluate  (Yk(N,a)}  up  to  k  =  J-l; 

from  equation  (C-9),  evaluate  {yj^N-I,!)}  up  to  k  =  J; 

from  equation  (C-ll),  evaluate  moments  {//y(k)}  up  to  k  = 
from  equation  (C-12),  evaluate  cumulants  (Xv(k)}  up  to  k 


-12) 

-13) 

9)) 

-14) 


J; 

=  J. 


C-3/(C-4  blank) 


APPENDIX  D  -  ROCS  FOR  KM  =  4,  RANDOM  GAUSSIAN  SIGNAL 

This  appendix  contains  the  ROCs  for  or-ing  with  pre-  and 
post-averaging,  when  the  time-bandwidth  product  KM  is  fixed  at  4; 
the  possible  combinations  (from  table  1)  are  repeated  here: 

KM  N 

4  1  1,2,4,8,16,32 

2  2 

1  4 


For  N  =  1,  only  the  product  KM  matters;  the  first  plot  in  this 
appendix  covers  this  special  case,  under  the  labeling  K  =  4, 

N  =  1,  M  =  1.  The  other  5  values  of  N,  along  with  the  3  possible 
combinations  of  K  and  M,  yield  15  additional  ROCs,  for  a  total  of 
16  ROCs  in  this  appendix. 


D-l 
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ROCs  for  K  =  4,  N  =  16,  M  =  1  Figure  D-6.  ROCs 
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ROCs  for  K  =  2,  N  =  8,  M  =  2  Figure  D-10.  ROCs 
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ROCa  for  K  =  1,  N  =  4,  M  =  4  Figure  D-14.  ROCs  for  K 


Figure  D-15.  ROCs  for  K  =  1,  N  =  16,  M  =  4  Figure  D-16.  ROCs  for  K 


APPENDIX  E  -  ROCS  FOR  KM  =  16,  RANDOM  GAUSSIAN  SIGNAL 


This  appendix  contains  the  ROCs  for  or-ing  with  pre-  and 
post-averaging  when  the  time-bandwidth  product  KM  is  fixed  at  16; 
the  possible  combinations  (from  table  1)  are  repeated  here: 

KM  N 

16  1  1,2,4,8,16,32 

8  2 

4  4 

2  8 

1  16 

For  N  =  1,  only  the  product  KM  matters;  the  first  plot  in  this 
appendix  covers  this  special  case,  under  the  labeling  K  =  16, 

N  »  1,  M  =  1.  The  other  5  values  of  N,  along  with  the  5  possible 
combinations  of  K  and  M,  yield  25  additional  ROCs,  for  a  total  of 
26  ROCs  in  this  appendix. 
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Figure  E-7.  BOCa  far  K  =  8,  N  =  2,  M  =  2  Figure  E-8.  ROCs  far  K 
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JROCb  far  K  =  8,  N  =  8,  M  =  2  Figure  E-10.  ROCs  far  K 
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Figure  B-27#  ROCs  for  K  ~  2,  N  ®  2,  M  —  8  Figure  E—18.  ROCs  for  K 
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Figure  E-23.  ROCs  for  K  =  1,  N  =  4,  M  *  16  Figure  E-24.  ROCs  far  K 


for  K  =  1,  N  =  16,  M  =  16  Figure  E-26.  ROCs  for  K 


APPENDIX  F  -  ROCS  FOR  KM  =  64,  RANDOM  GAUSSIAN  SIGNAL 


This  appendix  contains  the  ROCs  for  or-ing  with  pre-  and 
post— averaging  when  the  time— bandwidth  product  KM  is  fixed  at  64; 
the  possible  combinations  (from  table  1)  are  repeated  here: 


K 


M 


N 


64 

1 

32 

2 

16 

4 

8 

8 

4 

16 

2 

32 

1 

64 

1,2,4,8,16,32 


For  N  =  1,  only  the  product  KM  matters;  the  first  plot  in  this 
appendix  covers  this  special  case,  under  the  labeling  K  =  64, 

N  =  1,  M  =  1.  The  other  5  values  of  N,  along  with  the  7  possible 
combinations  of  K  and  M,  yield  35  additional  ROCs,  for  a  total  of 
36  ROCs  in  this  appendix. 


ROCs  for  K  =  64,  N  =  1,  M  =  1  Figure  F-2 .  ROCs 
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ROCs  for  K  =  32,  N  =  8,  M  =  2  Figure  F-10.  ROCs  for  K 
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Figure  F-13.  ROCb  for  K  =  16,  N  =  4,  M  =  4  Figuxo  F-14.  BOCs  for  K 


ROCs  for  K  =  16,  N  =  16,  M  b  4  Figure  F-16.  ROCs  for  K 
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Figure  F-17.  ROCs  far  K 
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ROCs  for  K  =  8,  N  =  32,  M  =  8  Figure  F-22.  ROCs  for  K 
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Figure  F-25.  ROCa  for  K  =  4,  N  =  16,  M  =  16  Figure  F-26.  BOCs  for  K 
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Figure  F-33 .  ROCs  for  K  =  1,  N  =  4,  M  =  64  Figure  F-34.  FOCs  far  K 


APPENDIX  6  -  ROCS  FOR  KM  =  256,  RANDOM  GAUSSIAN  SIGNAL 


This  appendix  contains  the  ROCs  for  or-ing  with  pre-  and 
post-averaging  when  the  time-bandwidth  product  KM  is  fixed  at 
256;  the  possible  combinations  (from  table  1)  are  repeated  here: 


K  M 


256 

1 

128 

2 

64 

4 

32 

8 

16 

16 

8 

32 

4 

64 

2 

128 

1 

256 

N 

1,2,4,8,16,32 


For  N  =  1,  only  the  product  KM  matters;  the  first  plot  in  this 
appendix  covers  this  special  case,  under  the  labeling  K  *  256, 

N  =  1,  M  =  1.  The  other  5  values  of  N,  along  with  the  9  possible 
combinations  of  K  and  M,  yield  45  additional  ROCs,  for  a  total  of 
46  ROCs  in  this  appendix. 
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Figure  G-l.  ROCs  for  K  =  256 ,  N  =  1,  M  =  1  Figure  G-2.  ROCs  for  K  —  256, 


Figure  6-5.  ROCa  for  K  =  256,  N  =  16,  M  =  1  Figure  G-6 .  ROCs  for  K  =  256 \ 
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Figure  6-9.  ROCa  for  K  =  128,  N  =  8,  M  =  2  Figure  6-10.  ROCs  for  K  =  128, 
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ROCs  for  K  =  64,  N  =  4,  M  =  4  Figure  G-14.  ROCs 


ROCs  for  K  =  64,  N  *  16,  M  =  4  Figure  G-16.  ROCs  for  K 
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Figure  G-23.  ROCs  for  K  =  16,  N  =  4,  M  =  16  Figure  G-24.  ROCs  for 
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Figure  G-37.  ROCs  for  K  =  2,  N  =  2,  M  =  128  Figure  G-38.  ROCs  for  K 
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G-22 


ROCs  for  K  =  1,  N  =  2,  M  =  256 
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Figure  G-45.  ROCs  for  K  =  1,  N  =  16,  M  =  256  Figure  G-46.  ROCs  for  K  =  1,  N  =  32,  M  =  256 


APPENDIX  H  -  STATISTICAL  DEPENDENCE  OF  0R-IN6 
OUTPUT  V  ( t )  AT  SEPARATED  TIME  INSTANTS 

The  initial  part  of  this  appendix  gives  the  results  of  a 
short  simulation  study  of  the  covariance  of  the  or-ing  output 
v(t)  at  time  separations  less  than  K,  while  the  second  part 
details  an  analytic  study  of  the  conditional  probability  of  the 
or-ing  output  v(t)  at  two  separated  time  instants. 

COVARIANCE  OF  OR-ING  OUTPUT 

The  particular  numerical  case  investigated  here  is  K  «■  4,  for 
which  the  pre-averager  outputs  are 

y  (t)  =  x  (t)  +  x  (t-1)  +  x  (t-2)  +  x  (t-3)  for  1  <  n  <  N,  (H-l) 
•*  n  n  n  n  n 

and  the  or-ing  output  is 

v(t)  =max{y1(t), - rYN(t)}  .  (H-2) 

The  results  of  a  simulation  for  hypothesis  Hq,  noise-only,  using 
50,000  independent  trials,  and  using  either  exponential,  flat,  or 
Gaussian  RVs  for  inputs  {xn(t)},  are  listed  in  table  H-l.  For 
special  case  N  =  1,  no  or-ing,  there  follows  v(t)  =  y^(t),  which 
is  a  sum  of  K  independent  inputs;  then,  v(t)  has  a  simple 
triangular  covariance  function,  going  linearly  to  zero  at  delay 
I  =  K.  This  effect  is  confirmed  by  the  N  =  1  simulation  results 
for  the  covariance  coefficient  in  table  H-l. 


H-l 


For  larger  N,  where  there  is  severe  competition  to  get 
through  the  "greatest-of "  device,  there  is  more  randomness 
between  adjacent  time  values  of  v(t).  This  causes  the  covariance 
coefficients  to  decrease  more  rapidly  with  the  time  separation  I. 
Then,  it  may  be  more  costly  to  increase  J  above  K/2  (roughly),  at 
least  for  large  N. 

Under  hypothesis  H^,  as  the  input  SNR  increases,  the  one 
signal-bearing  channel  tends  toward  dominance  of  the  greatest-of 
device;  therefore,  the  covariance  of  v(t)  should  tend  back 
towards  the  triangular  shape  as  the  input  SNR  increases. 

However,  since  the  processor  does  not  know  whether  HQ  or  ^  is 
prevalent,  it  may  be  operating  on  noise-only  samples;  this 
suggests  keeping  skip  factor  J  smaller  than  K  in  order  to  lessen 
losses. 


Table  H-l.  Estimated  Covariance  Coefficients  for  K  =  4 


I 

0 

1 

2 

3 

4 

N 

exponential 

1 

1 

.750 

.499 

.250 

.002 

2 

1 

.719 

.466 

.231 

.014 

4 

1 

.680 

.418 

.195 

.007 

8 

1 

.643 

.375 

.166 

.012 

1 

1 

.750 

.501 

.252 

.003 

flat 

8 

1 

.536 

.266 

.097 

-.005 

1 

1 

.749 

.499 

.248 

-.002 

Gaussian 

8 

1 

.553 

.278 

.104 

-.003 
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CONDITIONAL  PROBABILITY  OF  OR-ING  OUTPUT 


The  following  analysis  makes  no  assumptions  about  the 
statistics  of  the  inputs  {xn(t)}  except  independence.  However, 
the  particular  numerical  example  considered  will  presume  zero- 
mean  Gaussian  noise;  although  not  done  here,  the  extension  to 
signal-present  is  obvious.  The  parameters  M  and  J  are  not 
relevant  here,  because  post— averaging  is  not  yet  of  concern  at 
the  or-ing  output  v(t). 

The  pre-averager  outputs  in  the  n-th  channel  at  two  separated 
time  instants  K  and  K+I  are 

yn<K>  -  V1*  +  +  VK)  • 

V  (K+I)  »  X  (1+1)  +  •••  +  X  (K+I)  .  (H-3) 

n  n  n 

The  integer  variable  I  lies  in  the  range  0  <  I  <  K  for  this 
analysis.  Therefore,  the  covariance  coefficient  between 
pre-averager  outputs  yn(K)  an<*  Yn(K+I)  =  1  -  I/K. 

The  corresponding  or-ing  outputs  are 
v(K)  =  maxty^K),  ...  ,  YN(K)}  = 

=  max{x1(l)  +  •••  +  x^K),  ...  ,  xN(l)  +  •••  +  xN(K)}  (H-4) 

and 
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v(  K+l )  =  max  {y^  K+I),  ...  ,  yN(K+I)}  - 
=  maxfx^l+I)  +•••+  x^K+I),  ...  ,  xNU+I)  +•••+  xN(K+I)}  .  (H-5) 

At  this  point,  define  three  disjoint  subsequences 
An  -  xn(1)  +  •”  +  xn(I>  ' 

B  -  X  (1+1)  +  • • •  +  X  (K)  , 
n  n  n 

c  -  x  (K+l)  +  •••  +  x  ( K+I )  .  (H-6) 

n  n  n 

Since  the  inputs  are  independent  and  the  sequences  are  disjoint, 

every  RV  {Ar  Bx,  C1#  ...  ,  bn'  CN*  is  indePendent  of  everY 
other  one. 

Let  p„  and  c„  be  the  PDF  and  CDF,  respectively,  of  pre- 
averager  output  yR(t)  when  it  is  composed  of  K  terms,  as  in 
equation  (H-3).  Then,  the  joint  cumulative  probability  of  the 
two  or-ing  outputs  (H-4)  and  (H-5)  is 

Pr{v(K)  <  a,  v(  K+I )  <  (3}  - 

=  Pr{A^+B^<a,  .  .  •  ,  A^+Bj|<a,  B^+C^<S,  •  •  •  ,  — 

=  [Pr{A1+B1<a,  Bj+Cj^S}]”  =  I  Jdu  pK_x(u)  Cj(o-u)  Cj(p-u)]N.  (H-7) 

At  the  same  time,  the  ordinary  first-order  CDF  is 

Pr{v(K)<a}  =  [cR(a)]N  .  (H-8) 

Therefore,  the  conditional  probability  of  interest  is 
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Pr {v( K+I ) <0 | v(K) <a}  - 


rx  du  pR_I(u)  CjCa-u)  CjO-u)"! 


cR(a) 


N 


N 


=  [Q(a,P)] 

( H-9 ) 


This  physically-meaningful  measure  of  dependence  needs  only  one 
integral  evaluation  once  PDF  pR  and  CDF  cR  are  available.  (By 
contrast,  the  covariance  coefficient  of  v(t)  requires  much  more 
numerical  effort.)  Notice  that  this  conditional  probability  goes 
to  zero  as  N  +  ®,  that  is,  as  the  amount  of  or-ing  increases, 
regardless  of  the  value  of  separation  I  >  0.  Thus,  increased 
competition  at  the  or-ing  input  leads  to  more  independent  or-ing 
outputs,  even  for  1=1,  an  adjacent  time  sample. 


As  checks  on  this  result,  for  1  =  0,  then  AR  =  0,  CR  =  0, 
Pg(u)  =  S(u),  cQ(u)  =  U(u) ,  and  there  follows 


Q(a,e) 


c  (min{a,{3} ) 


1 


[cK(P)/cK(a) 


\ 

for  0  >  a 
for  0  <  a 

/ 


( H-10 ) 


On  the  other  hand,  for  I  =  K,  there  follows 

c_ ( a)  c_0) 

p^tu)  =  Pq ( u )  =  8(u)  ,  Q(a,e)  =  R  c  (ay -  =  cR(p)  .  (H-ll) 

K 


Both  of  these  special  cases  are  obviously  correct. 


As  a  numerical  example,  consider  zero-mean  unit-variance 
Gaussian  noises  {xn(t)}.  Then, 
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Pj(u)  = 


1 


(  2  ill ) 15 


ci <u ’  *  *(tt)  ' 


-  Jrwm  J  dt  ♦<t)  *(~  </iI>'8t)  -)  •  (H-12) 


A  couple  of  numerical  examples  of  conditional  probability  (H  9) 
follow,  which  confirm  the  increased  statistical  dependence  as 
separation  I  tends  to  zero. 


a-1  f  (3=1,  K=4 : 


I 

1 

2 

3 

4 

QU,3) 

.854 

.790 

.738 

.691 

ot —  2  ,  8 — 2 , 

I  I 

K=10 : 

1  2 

3 

4 

5 

6 

7 

8 

9 

10 

Q(a,(3)  | 

.920  .887 

.861 

.838 

.818 

.800 

.783 

.767 

.751 

.736 
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APPENDIX  I  -  MATLAB  PROGRAM  FOR  EVALUATION 
OF  ROCS  FOR  RANDOM  GAUSSIAN  SIGNAL 

The  parameters  of  interest  are  entered  at  the  top  of  the 
program  listed  here.  The  sampling  increment  du  is  applied 
directly  to  the  PDF  of  v(t)  and  controls  the  spacing  of  the 
aliasing  lobes  in  the  calculated  CF  of  v(t).  If  the  lobes  are 
too  close,  yielding  significant  aliasing  in  the  CF  domain,  the 
value  of  du  must  be  decreased.  The  range  of  input  SNR  values  to 
be  investigated  is  governed  by  rmin  and  rinc. 

For  each  input  SNR  value  (isnr),  the  magnitude  of  the  CF  of 
v(t)  is  first  plotted,  so  that  the  aliasing  lobes  can  be  viewed; 
if  the  lobes  are  too  close,  du  must  be  decreased  and  the  CF 
rerun.  When  satisfactory  CF  lobe  spacing  has  been  achieved, 
calculations  continue  with  the  corresponding  EDF  of  w(t),  which 
is  plotted.  If  aliasing  is  a  problem  in  this  EDF  domain,  the 
size  n  of  the  FFT  must  be  increased.  Finally,  when  satisfactory 
values  of  du  and  n  have  been  realized  for  all  values  of  SNR  of 
interest,  the  ROCs,  namely,  Pd  versus  Pf,  are  plotted  for  all  the 
input  SNR  values  considered. 

To  find  the  required  input  SNR  necessary  to  achieve  a 
specified  and  Pd,  these  probability  values  must  then  be 
entered,  along  with  the  number  of  the  ROC  that  is  closest  to  the 
operating  point  of  interest.  Linear  interpolation  is  conducted, 
and  the  required  input  SNR  is  printed  out. 
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clear,  elf  %  NM2  TR  11,150 

b=0;  %  Additive  bias  to  w 

du=.l;  %  Sartpling  incranent  in  u 

n=2Al 2;  %  EFT  size 

kc=4;  %  K,  amount  of  pre-averaging 

nc=l;  %  N,  amount  of  or-ing 

irc=l;  %  M,  amount  of  post-averaging 

rmin=-3;  %  Mmimum  SNR  (dB) 

rinc=l;  %  Incranent  in  SNR  (dB) 

num=17;  %  Nunber  of  ROCs 

disp(  'b  du  n;  kc  nc  me: ' ) 

disp  ( [b  du  n;  kc  nc  me] ) 

xg=[le-6  le-5  le-4  .  001  .002  .  005  .  01  .02... 

.05  .1  .2  .3  .4  .5]; 

yg=[le-6  le-5  le-4  .001  .002  .005  .01  .02... 
.05  .1  .2  .3  .4  .5  .6  .7  .8  .9  .95  .98  .99]; 
jg=phiinv(xg)  ; 
yg=phiinv(yg)  ; 

[Xg,Yg]  =meshgrid  (xg,yg)  ; 

Rfegcf=zeros  (n,  1)  ; 

EDF=zeros(n,l) ; 

Bd=zeros  (n,nurm-l)  ; 

nl=n-l;  n2=n/2;  kcl=kc-l;  ncl=nc-l;  n3=n2+l; 

pn=pi/n;  d2=du/2; 

arg=pn.*(l:n2)  '  ; 

sinc=sin  (arg) .  /arg; 

dispC  isnr  edfO  -  1') 

for  isnr=0:num 

rhodb=rrnin+rinc*  (isnr-1)  ; 
rho=10A  ( .  l*rhodb)  ; 
if  (isnr=0)  rho=0;  end 
as=l/  (1+rho)  ; 

X=zeros  (n,  1)  +i*zeros  (n,  1)  ; 
meanv=0;  cdfk=0;  k=-l; 
vMle(cdfk<.5  |  area>le-20) 
k=k+l;  uk=du*k;  u2=uk+d2; 
au=as*u2;  p0=exp(-u2) ;  e0=p0; 
pl=exp  (-an) ;  el=pl; 
for  j=l:kcl 

p0=p0*u2/j;  e0=e0+p0; 
pl=pl*au/j;  el=el+pl; 
aid 

cO=max  (1-eO ,  0) ;  cl=max  ( 1-el ,  0 )  ; 
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cdfo=cdfk;  cx3fk=cl*cOAncl; 

area=odfk-cdfo; 

j=mod(k,n); 

X  ( j+1)  =X(  j+1)  +area; 
iteanv=meanv+area*uk;  %  uk,  not  u2 

aid 

X=fftgreen(X)  ; 

Magcf =loglO  (X.  *conj  (X)  +le-50)  * .  5; 
plot  (Magcf) 

axis([l  n+1  -16  0]);  grid  on 
pause (1) 

cbd=2*pi/  (n*du)  ; 

X(2:n3)=conj  (X(2:n3) )  ./sinc(l:n2)  ; 
X(n2+2:n)=0; 

X=X./'irc; 

X(l:n3)=X(l:n3)  .*exp(i*b*dxi* (0:n2)  •),- 
neanw=3iieanv^itc+b; 

X(1)=0; 

X(2:n3)=X(2:n3)  ./(l:n2)  '; 

X=fftgreen(X)  ; 
a=.5-nneanw/  (n*du)  ; 
edfO=a+imag(X(l) )  /pi; 
disp([isnr  edfO-1]) 
k=(l:n) 

edf=a- (k-1)  ./n+imag(X(k) )  ./pi; 

X=edf; 

EDF=loglO  (abs  (edf )  +le-30)  ; 
plot (EDF) 

axis  ( [1  n+1  -18  0] ) ;  grid  on 
paused) 
edf=real(X(k) )  ; 
edf=min(edf ,  l-le-12 )  ; 
edf =max  (edf ,  le-12 )  ; 

Bd( : , isnr+1)  =phiinv(edf) ; 
end 
beep 
pause 
elf 

hold  cn 

set(gcf, 'RaperPosition' ,  [.25  .25  8  10.5]) 
plot(xg,Yg,  'k‘ ) 
plot  (Xg,yg,  *k* ) 

plot  (xg,xg,  'k')  %  zero  SNR  ROC 


plot (Bd( : ,  1)  [2:num+l] ) ,  'k') 

axis([xg(l)  xg(14)  yg(l)  yg(21)]) 
axis  off 
vMle  1 

T=input  ( 'pf  pd  isnr:  ' )  ; 
if  T(1)=0  break;  end 
pft=phiinv(T(l) )  ; 
pdt=piIiiinv'(T(2) )  ; 
isnr=T(3) ; 
for  j=0:nl 

if  (Bd(j+l,l)<pft)  break;  end 
end 

xl=Ed(j+l,l)  ; 
x2=Ed(j,l); 
as=Ed  ( j +1 ,  isnr+1 )  ; 
bs=Pd(j+l, isnr+2) ; 
cs=Bi(j,  isnr+1)  ; 
ss=  (cs-as)  /  (x2-xl)  ; 
al=pdt-ss*pft; 
fs=  (al+ss*xl-as)  /  (bs-as)  ; 
rc=rmin+rinc*  (isnr-l+fs)  ; 
disp(  [T  fs  rc] ) 
ad 

%  function  y=phiinv(x)  for  0  <  x  <  1 
%  y=1.414213562373095*erfinv(2*x-l)  ; 
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